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ABSTRACT 


A rapid  and  accurate  approximation  to  the  Lorentz  line  molecular 
absorption  coefficient  has  been  developed  which  significantly  reduces 
the  time  of  line-by-line  transmittance  calculations.  By  separating 
computations  required  for  each  spectral  line,  atmospheric  level,  and 
absorbing  gas  under  consideration,  arithmetical  operations  are 
minimized . 

Time  tests  were  conducted  while  computing  the  monochromatic 

absorption  coefficient  for  five  frequency  channels  in  the  4.3  hm' 

band.  The  approximation  was  compared  to  the  evaluation  of  the 

Lorentz  line  shape  (standard  manner)  for  a 33  level  atmosphere  and 

an  average  of  1000  absorption  lines  per  region^  Water  vapor,  carbon 

dioxide,  nitrous  oxide,  and  carbon  momoxide  molecules  »:ere  considered. 

The  approximation  was,  in  all  cases,  at  least  9 times  faster. 

> Absorption  coefficient  values  computed  by  this  approximation 

agreed  to  those  computed  in  the  standard  fashion  to  at  least  three 

decimal  places  regardless  of  molecule  type,  channel,  or  atmospheric 

-4  ' 

level.  Average  percent- relative  errors  were  typically:  10  for 

.1)00'  --- 

H,0;  10'4  for  C0~:  10"3  for  N,0;  10-8  for  CO. 

L ,atM  1 ■ .ocoosOCl 

For  transmittance  calculations  over  a Land,  time  savings  can  be 
predicted  by  analyzing  the  required  number  of  arithmetical  operations. 
For  1,000  spectral  lines  over  a bandwidth  of  25  cm and  a step  in- 
crement of  0.02  cm  a reduction  in  computation  time  of  450Z  is 


anticipated. 
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CHAPTER  I; 


INTRODUCTION 

One  of  the  most  widely  used  methods  of  calculating  atmospheric 
abaorptim  Is  the  line-by-line,  direct  integration  method.  With  the 
aid  of  modern  high-speed  computers  and  the  availability  of  accurate 
line  parameters,  an  ever  increasing  confidence  is  being  placed  in 
this  technique.  The  adversion  to  it,  however,  is  that  computation 
time  and,  consequently  cost,  tend  to  excalate  beyond  feasible  limits. 
Even  though  larger  and  faster  computers  are  constantly  being  developed, 
the  need  still  exists  for  a fast  and  accurate  method  of  calculating 
the  molecular  absorption  coefficient. 

About  a year  ago,  Dr.  J.  H.  Pierluissi  122]  and  Dr.  G.  A.  Gibson 
110] , under  a research  contract  from  the  Atmospheric  Sciences  Labora- 
tory, WSMR,  began  investigating  a method  of  reducing  the  number  of 
required  calculations  to  obtain  a molecular  absorption  coefficient. 

This  work  will  elaborate  on  their  proposed  approximation  to  the  Lorentz 
line  shape  which  promises  to  significantly  reduce  the  time  and  cost 
of  line-by-line  calculations.  Because  application  of  this  approxima- 
tion to  the  band  problem  is  not  final,  only  the  monochromatic  case  will 
be  fully  developed.  However,  a discussion  of  polychromatic  radiation 
la  Included.  Evidence  of  reduction  in  the  number  of  required  calcula- 
tions in  the  standard  Lorentz  line  shape  is  offered  in  Chapters  5 and 
6. 

for  the  sake  of  completeness , a chapter  has  been  devoted  to  the 
general  theory  of  molecular  absorption.  Additionally,  some  space  has 


1 


been  allocated  to  the  "AFCRL  Atmospheric  Absorption  Line  Parameters 
Compilation"  data  tape  jMcClatchey^J , which  was  used  in  all  calcula- 
tions and  to  the  units  generally  used  in  atmospheric  sciences  (Chapter 
4). 

All  computer  programs  presented  were  written  in  Fortran  IV  and 
designed  to  be  used  with  the  University  of  Texas  at  El  Paso's  IBM- 
360/65,  G-level  compiler.  Some  slight  modification  of  the  job-control 
cards  may  have  to  be  made  if  these  programs  are  to  be  used  on  another 
computer. 
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CHAPTER  II 

GENERAL  THEORY  OF  ATMOSPHERIC  TRANSMITTANCE  AND  ABSORPTION 

This  chapter  will  furnish  the  reader  with  a basic  review  of  the 
fundamentals  of  atmospheric  transmittance  of  radiant  energy.  The 
first  section  has  been  devoted  to  the  earth's  atmosphere  with  empha- 
sis placed  on  its  classification  methods  and  composition.  Section 
2.3  is  a discussion  of  molecular  energies  and  line  spectra  resulting 
from  state  transitions.  This  is  followed  by  an  explanation  of  the 
different  spectral  line  broadening  mechanisms.  Finally,  the  chapter 

concludes  with  the  formulation  of  the  Lambert-Beer  law  for  molecular 
4 

absorption  [Brugel  ] and  a discussion  of  the  principle  gases  involved. 
2.1  Earth's  Atmosphere 

There  are  many  methods  for  classifying  the  earth's  atmosphere, 
but  the  most  widely  accepted  method  is  based  on  temperature  varia- 
tions with  altitude.  One  such  distribution  is  presented  in  Fig.  2.1-1. 
A region  of  the  atmosphere  is  characterized  by  its  temperature  gradi- 
ent (either  positive  or  negative)  and  has  the  suffix  "-sphere."  As 
the  variation  in  temperature  begins  to  change  sign,  indicating  the 
beginning  of  a new  region,  the  suffix  is  changed  from  "-sphere"  to 
"-pause"  (i.e.  stratosphere  to  stratopause) . 

The  lowest  region,  0-15  km,  is  known  as  the  troposphere.  It  is 
characterized  by  a negative  vertical  temperature  gradient.  The  ceil- 
ing, tropopause,  may  vary  as  much  as  10  km  depending  on  geographical 

location  and  season.  Temperatures  typically  decrease  6°K  per  km 
3 

[ Banks  ] . 
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Temperature  in  the  "stratost  here,"  15-50  km,  increases  to  a 


maximum  of  about  270°K.  The  hearing  effect  is  due  to  the  absorption 

of  ultraviolet  radiation  by  ozone  (0^)  and  to  atmospheric  infrared 

3 

emissions  [Banks  ] . 

The  "mesosphere"  rests  above  50  km  and  reaches  to  a height  of 
about  100  km.  Here  the  temperature  reaches  a minimum  at  85+5  km. 

The  highest  region  shown  in  Fig.  2.1-1  is  the  "thermosphere." 

It  extends  beyond  100  km  and  has  a positive  temperature  gradient. 
Typically,  unless  specific  requirements  dictate  otherwise,  classifi- 
cation of  the  atmosphere  can  stop  with  the  "thermosphere." 

While  temperatures  normally  vary  up  or  down  with  altitude, 

atmospheric  pressure  always  decreases.  Variations  in  pressure  are 

28 

related  to  those  in  altitude  by  the  hydrostatic  equation  [Whitten  ] 

dP  = -g(z)  p(z)  dz,  2.1-1 

where  P is  the  atmospheric  pressure,  g(z)  is  the  acceleration  due  to 
gravity  at  altitude  z,  and  p(z)  is  the  atmospheric  density  at  alti- 
tude z.  Fig.  2.1-2  illustrates  the  changes  in  pressure  with  respect 
to  altitude. 

In  terms  of  gas  concentrations,  the  atmosphere  is  composed  of 
two  groups  of  molecules.  Those  which  are  mixed  uniformly  in  the 
atmosphere  are  listed  in  Table  2.1-1.  Those  which  have  a variable 
mixing  ratio  are  water  vapor  0^0)  and  ozone  (0^) . The  densities  of 
these  molecules  (H^O,  0^)  in  the  atmosphere  are  highly  dependent  on 
geographical  location,  season,  and  even  time  of  day. 
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ATMOSPHERE 

CHEMICAL 

PERCENT 

CONSTITUENT 

FORMULA 

BY  VOLUME 

ABSORBER 

Nitrogen* 

*2 

78.084 

No 

Oxygen* 

°2 

20.946 

No** 

Argon* 

A 

0.934 

No 

Carbon  Dioxide* 

C02 

0.032 

Yes 

Neon 

Ne 

1.818  x 10"3 

No 

Helium 

He 

5.24  x 10~4 

No 

-4 

Methane 

ch4 

2.0  x 10 

Yes 

-4 

Krypton 

Kr 

1.14  x 10 

No 

-5 

Nitrous  Oxide 

n2° 

5.0  x 10 

Yes 

-5 

Hydrogen 

H2 

5.0  x 10 

No 

-6 

Xenon 

Xe 

9.0  x 10  0 

No 

*Major  atmospheric  constituent 

**McClatchey , et.  al.,  has  included  490  lines  of  O2  in  AFCRL's 
data  compilation  [17]. 


TABLE  2.1-1  COMPOSITION  OF  THE  EARTH'S  ATMOSPHERE  (DRY) 

12 

[Adapted  from  Hudson  ] 


Another  method  for  categorizing  the  contents  of  the  atmosphere 
is  by  classifying  each  component  as  a major  or  minor  constituent. 
The  major  constituents,  Nitrogen  (N2) , Oxygen  (C>2),  Argon  (A),  and 
Carbon  Dioxide  (CO^)  , make  up  over  99.9  percent  of  the  atmosphere. 


2.2  Molecular  Energy  and  Spectra 

Because  molecular  and  atomic  spectra  both  exhibit  discrete 
energy  transitions,  the  fundamental  atomic  assumptions  made  by  Niels 
Bohr  in  1913  [Fowles^]  can  be  used  to  illustrate  the  quantization  of 
molecular  energies.  In  a molecular  sense,  these  assumptions  tell  us 
that  molecules  can  possess  only  certain  discrete  energy  levels.  When 
a transition  from  one  state  to  another  occurs,  the  frequency  at  which 
the  transition  takes  place  is  governed  by  Plank's  relation.  This 
frequency  (v)  is  given  by 


where  AE  is  the  energy  difference  between  the  two  states  involved  and 
h is  Plank's  constant. 

Fig.  2.2-1  illustrates  two  possible  transitions  between  energy 
states  and  E2<  A rise  in  energy  from  E^  to  indicates  that  the 
molecule  has  "absorbed"  incident  radiation.  Similarly,  a lowering  in 
energy  may  be  described  as  an  "emission." 

In  general,  the  total  energy  of  a molecule  is  a function  of  four 
factors  and  their  combinations.  The  first  is  the  energy  due  to  trans- 
lational motion  of  the  molecule  which  is  dependent  on  velocity.  For 
the  purposes  of  molecular  absorption  analysis  the  translational 
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a.  Molecular  energy  is  raised  from  lower  E^  state  to  E2  by  absorbing 


radiant  energy. 

b.  Molecular  energy  is  lowered  from  higher  Ej  state  to  E^  by  emit- 
ting radiant  energy. 


Figure  2.2-1  STATE  ENERGY  TRANSITIONS 
[Adapted  from  Fowles^] 
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energy  can  be  ignored  iGoody^j  . 

Hence,  this  work  only  considers  three  types  of  energy  which 
contribute  to  total  energy.  They  are  the  electronic  energy  (E^), 
the  rotational  energy  (Erot)>  and  the  energy  due  to  the  vibrational 
modes  of  motion  (E^^) . Therefore, 


tot 


rot 


+ E 


vib 


+ E , + E 


el  comb  * 


2.2-2 


where  E , represents  the  interaction  of  two  or  more  energy  types, 
comb 

The  largest  contributor  to  total  energy  is  that  of  the  electronic 

type  which  is  typically  a few  electron-volts.  Vibrational  energies 

are  a few  tenths  or  hundredths  of  an  electron-volt,  and  the  smallest 

is  due  to  rotation  which  is  on  the  order  of  hundredths,  thousandths 

29 

or  ten- thousandths  of  an  electron-volt  [Zeuv  ] . 

Each  of  the  three  energy  types  are  quantized  in  a different  way 
and  hence,  by  Eq.  2.2-1,  may  absorb  or  emit  at  different  frequencies. 
Transitions  that  involve  rotational  energy  levels  only  take  place  in 
the  far  infrared  or  the  microwave  region  of  the  spectrum.  It  is 
known  as  the  "pure  rotation  spectrum."  A "rotational-vibrational 
spectrum"  is  also  possible  and  can  be  found  in  the  near  infrared. 

The  most  energetic  transitions  are  due  to  electronic  energy  changes. 
The  "electronic  spectra"  for  molecules  occur  in  the  visible  and 
ultraviolet  region.  See  Fig.  2.2-2. 

The  scope  of  this  document  will  not  allow  an  elaborate  expla- 
nation of  Quantum  Theory.  It  will  suffice  to  state,  however,  that 
the  frequencies  at  which  a molecule  can  absorb  or  emit  radiant 
energy  can  be  predicted  from  this  theory.  This  is  done  by  assigning 


i 
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"quantum  numbers"  to  each  of  the  molecule's  energy  types.  These 
numbers  then  describe  the  allowable  energy  transitions  or  allowable 
frequencies  for  emission  or  absorption.  The  assignment  of  "quantum 
numbers"  is  based  on  the  individual  molecule's  geometry,  modes  of 
motion,  and  number  of  atoms. 

The  absorption  or  emission  frequencies  which  are  predicted  by 
this  theory  are  known  as  spectral  lines.  The  strength,  line  width, 
and  energy  associated  with  each  spectral  line  for  a given  molecule 
is  the  subject  of  much  research  and  study.  A widely  used  compila- 
tion of  these  parameters  by  R.  A.  McClatchey,  et.  al.,  [17]  is  dis- 
cussed in  chapter  4. 

2.3  Line  Broadening  Mechanisms 

Unfortunately,  when  a transition  between  energy  states  takes 
place  in  a molecule  the  corresponding  spectral  line  can  not , in  real- 
ity, be  considered  to  be  purely  monochromatic.  Depending  on  atmo- 
spheric conditions,  spectral  lines  exhibit  a broadened  contour  shape. 
The  broadening  can  be  attributed  to  the  effects  of  radiative  damping, 
doppler  shifts,  and  molecular  collisions. 

Natural  spectral-line  broadening  is  a result  from  the  fact  that 
energy  transitions  take  a finite  amount  of  time  to  occur.  Hence,  the 

vibrational  process  that  accompanies  the  transition  gives  rise  to  a 

29 

spectral  interval  [Zuev  ].  Excluding  all  other  disturbances  the 
natural  line  shape  is  given  by 


l 


12 


where  the  line  width  is  given  as 


Y 


N 


1 

2irt 


2.3-2 


vq  is  the  center  frequency  of  the  line,  and  t is  the  time  spent  by 

the  molecule  in  the  excited  state. 

The  form  of  Eq.  2.3-1  will  be  seen  again  in  the  Lorentz  line 

shape  and  in  this  sense  natural  broadening  is  important.  However,  the 

effects  are  normally  ignored  in  transmittance  calculations  because  of 

the  small  line  width.  Typically  is  on  the  order  of  10  ^ wave  num- 
29 

bers  [Zeuv  ] . 

In  the  regions  of  the  atmosphere  where  pressures  are  around  10  mb 
and  below,  the  line  broadening  effects  due  to  the  doppler  principle 
are  dominant  [Drayson^].  Assuming  a maxwellian  velocity  distribution 
among  the  molecules  and  ignoring  all  other  effects,  the  Doppler  line 
shape  is  predicted  by 

2 

K(v)  * kQ  e"y  , 2.3-3 

where 


y 


Iv-vQ]  Iln2] 


1/2 


2.3-5 


and  aD  is  the  Doppler  half-width. 


r 

I * 

The  line  strength,  S,  is  a function  of  temperature  and  is  expressed 

' s-s°(r)  (*vBtp  ir-  (f2-1) 

where  b is  a gas  dependent  constant,  is  the  vibrational  partition 
function  ratio,  E"  is  the  rest  energy,  and  k is  Boltzmann's  constant. 

It  should  be  noted  here  that,  in  most  texts,  Eq.  2.3-6  is  given 
without  the  ratio  of  partition  functions.  In  this  report,  the  ratio 
will  always  be  included  so  that  calculations  involving  S will  agree 
more  closely  with  the  line  parameter  data  that  was  used  [McClatchey^] . 

As  atmospheric  pressure  increases,  molecular  collisions  become 
much  more  frequent.  Hence,  it  is  natural  to  suspect  that  the  line 


broadening  mechanism  in  such  regions  of  the  atmosphere  is  heavily  de- 
pendent on  pressure  (molecular  collisions) . The  most  widely  accepted 
formulation  of  the  pressure  broadened  line  shape  is  attributed  to 

H.  A.  Lorentz  (1906)  although  works  in  this  same  area  were  done  earlier 

28 

by  A.  A.  Michelson  (1895)  and  later  by  W.  Lenz  (1924)  [Goody  ].  This 
important  line  shape  is  discussed  below. 

As  mentioned  in  the  beginning  of  this  section,  the  spectral  line 
shape  describes  how  a theoretically  single  spectral  line  is  actually 
broadened  by  various  mechanisms.  The  deviation  from  mono chromatism 
can  be  factored  into  the  product  of  a line  strength  and  a line  shape. 
This  line  shape  is  usually  normalized  so  that  its  integral  over  infi- 
nite limits  is  unity. 

Assuming  that  the  Pressure-broadened  line  shape  could  be  based 


14 


on  "strong"  encounters  between  molecules  and  a statistical  function 
describing  the  time  between  collisions,  Lorentz  developed  the  well- 
known  formulation 


- (i) 


2.3-7 


Coupling  the  line  shape  and  line  strength  (Eq.  2.3-6)  we  have 


k(v)  - f 


it  2 . .2  * 

aT  + (v-v  ) 

L o 


2.3-7a 


The  most  simplified  expression  for  the  Lorentz  half-width  used 
in  atmospheric  applications  is 


r p irT°i  1/2 

°L  " “L  L^JLt-J 


2.3-8 


where  a is  the  half-width  value  at  STP.  A contour  of  the  Lorentz 
o 

profile  is  shown  in  Fig.  2.3-1. 

When  transmittance  calculations  are  to  be  performed  over  a wide 
range  of  pressure  values,  a combined  Doppler-Lorentz  contour  is  some- 
times used  to  gain  better  accuracy.  The  convolution  of  these  two  line 
shapes  is  known  as  the  Voigt  profile  [La  Rocca^  ] . It  is  expressed  as 

2 

dt,  2.3-9 


k U 


k(v) 


-t 


2 2 
IT  + (y-t; 


where  kQ  and  y are  defined  in  Eqs.  2.3-4  and  2.3-5  respectively,  and 
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MAGNITUDE 


Figure  2.3-1  LORENTZ  LINE  SHAPE 
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« 


t is  a dummy  yariable.  U,  a function  of  the  Doppler  and  Lorentz  line- 
widths,  is 


U 


aL  Iln2]1/2 


2.3-10 


It  is  of  Interest  to  note  that  while  the  Voigt  profile  delivers 
greater  accuracy  over  wide  pressure  ranges  (i.e.  0-1013  mb  for  earth's 
atmosphere),  it  is  difficult  and  time  consuming  to  evaluate.  This  is 
due  primarily  to  the  integration  over  infinite  limits  and  the  vari- 
able characteristics  of  the  Doppler  half-width  [Goody1^] . An  interest- 
ing work  in  this  area  has  been  compiled  by  Van  Derwood  [26], 

2.4  Beer's  Law  of  Molecular  Absorption 

Attenuation  of  radiant  energy  in  the  atmosphere  can  be  attributed 
to  the  process  of  scattering  and  absorption.  The  latter,  because  of 
allowable  state  transitions,  is  the  primary  cause  of  such  abatements 
in  the  infrared  spectrum.  A quantitative  analysis  of  molecular  absorp- 
tion is  given  by  Beer's  law  (also  known  as  the  Lambert-Beer  law)  which 
is  outlined  below. 

Referring  to  Fig.  2.4-1,  let  AI  (z)  be  the  amount  of  radiation 

o 

which  is  absorbed  by  the  gas  contained  within  Az.  The  subscript  indi- 
cates monochromatic  radiation  at  frequency  vq.  The  gas  amount,  U(z), 
between  z^  and  z ^ can  be  computed  by 

U(z)  - p(z)Az  2.4-1 

if  the  gas  density  is  known  in  the  layer. 

i 
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Figure  2.4-1  LAMBERT-BEER  LAW  GEOMETRY 
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The  absorption  law  of  Beer  states  that  the  loss  in  intensity 
between  z^  and  is  proportional  to  the  path  length,  Az,  the  incom- 
ing radiation,  I (z^) , and  the  concentration  of  absorbing  gas,  p(z). 
o 

If  K (z)  represents  the  constant  of  proportionality,  then 
o 

AIv  (z)  - -Iv  C*l)Kv  U)p(z)Az  . 2.4-2 

o o o 

Letting  the  path  length  become  arbitrarily  sniall,  Eq.  2.4-2  can  be 
written  in  differential  form  as 


dl  (z)  = -I  (z)Kv  (z)p(z)dz, 
o V0  o 


2.4-3 


where  the  negative  sign  indicates  the  loss  or  absorption  of  radiant 

intensity.  Letting  t be  the  ratio  of  incoming  to  outgoing  intensi- 

o 

ties,  and  integrating  Eq.  2.4-3  between  and  z^  gives 


T 


V 

o 


Z2 

- EXP 

- 

K (z)p(z)  dz 

V 

0 

Z1 

2.4-4 


If  the  gas  density  and  constant  of  proportionality  are  known,  Eq.  2.4-4 
becomes 


t - exp 

v 

o 


-K  U 

v 

o 


2.4-4 


where  the  gas  amount,  U,  is  defined  by  Eq.  2.4-1. 

The  constant,  K , is  the  absorption  coefficient.  Sometimes 

Vo  « 
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t 


Is  referred  to  as  an  "extinction"  coefficient,  although  the  term  is 

normally  reserved  to  describe  a combination  of  attenuation  mechanisms. 

K is  formed  from  the  product  of  a given  line  shape  as  described  in 
o 

section  2.3  and  a line  strength  as  given  by  Eq.  2.3-6. 

Since  t measures  how  well  radiant  energy  passes  through  a 
o 

given  amount  of  absorbing  gas,  it  is  known  as  transmittance.  Trans- 
mittance ranges  from  values  of  1 (indicating  no  absorption)  to  0 (in- 
dicating total  absorption).  Additionally,  the  term  absorptance  applies 
to  a medium  abating  Incident  radiation  and  is  related  to  transmittance 
by  A ■ 1— t . Therefore,  using  Eq.  2.4-5,  absorptance  can  be  expressed 
in  terms  of  an  absorption  coefficient  as 


A = 1 - EXP 


-K  U 

v 

o 


2.4-6 


2.5  The  Absorbing  Gases 

As  noted  earlier,  most  molecular  absorption  takes  place  in  the 
Infrared  spectrum.  This  chapter  concludes  with  a brief  discussion 
of  the  atmospheric  constituents  which  are  mainly  responsible  for  this 
phenomenon . 

It  is  generally  accepted  that  seven  molecules  contribute  to 
absorption  in  the  infrared  and  visible  regions  of  the  atmosphere 
[2,  12,  16,  18  & 29].  They  ares  (1)  water  vapor  - H20;  (2)  car- 
bon dioxide  - C02:  (3)  ozone  - C>3;  (4)  nitrous  oxide  - N20;  (5)  car- 
bon monoxide  - CO;  (6)  methane  - CH^;  and  (7)  oxygen  - 02«  McClatchey, 
et.  al.,  considers  only  these  constitutents  in  AFCRI’s  "Line  Para- 
meter Compilation"  tape  which  is  the  source  data  for  all  calculations 
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used  in  this  document 


Of  the  seven  molecules,  ar'.y  CO ^ and  0^  are  major  atmospheric 
constituents.  Luckily,  oxygen  is  a very  minor  absorber.  The  other 
two  major  constituents,  N.,  and  A,  are  transparent  in  the  spectral 
region  of  interest.  Carbon  dioxide  is  the  only  absorber  which 
constitutes  any  significant  portion  of  the  atmosphere  (0.032%). 

Absorptance  calculations  are  most  difficult  for  H^O  and  0^  due 
to  their  variable  mixing  ratios,  Concentrations  are  dependent  on 
geographical  location  and  time  of  year.  Standard  density  profiles 

for  these  two  gases  are  available  for  typical  latitudes  and  seasons 

18  i 5 

[McClatchey  , LaRocca'  ] , 

The  other  molecules  are  considered  to  be  uniformly  mixed  in  the 
atmosphere.  Their  parts  per  million  by  volume  (PPM)  ratios  are  shown 
below. 


Constituent 

Ratio  in  PPM 

C02 

330.0 

n2o 

0.28 

CO 

0.075 

ch4 

1.60 

°2 

2,10  x 105 

[McC latch-  / ‘ | 

It  is  possible  i:o  ru'dJyze  each  molecular  gas  to  establish  the 
lines  of  absorption  ovei  i range  of  frequencies.  For  the  portion  of 
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I 

the  electromagnetic  spectrum  from  0-15  pm,  Fig.  2.5-1  illustrates 
the  percent  absorption  for  each  molecule  of  interest.  Each  profile 
is  characterized  by  discrete  bands  consisting  of  a distribution  of  a 
large  number  of  overlapping  lines  of  various  strengths.  The  overall 
contours  as  shown  in  Fig.  2.5-1  are  dependent  on  each  molecule's 
distribution  of  lines,  each  line's  half --width  and  strength,  and 
proximity  to  another  line. 


CHAPTER  III 


A PROPOSED  APPROXIMATION  TO  THE  LORENTZ  LINE  SHAPE 

3 . 1 Preface 

The  use  of  transmittance  calculations  based  on  remote  radiation 
measurements  to  predict  temperature  and  water  vapor  profiles  was 
first  suggested  in  the  late  1950's  by  Dr.  Lewis  D.  Kaplan  [15].  In 
order  to  retrieve  such  information  from  satellite  data  in  a timely 
fashion,  the  solution  of  the  transmittance  equation  (Eq.  2.4-4)  must 
be  accomplished  a3  accurately  and  as  rapidly  as  possible.  A proposed 
method  of  calculating  an  absorption  coefficient  based  on  the  Lorentz 
line  shape  is  outlined  in  this  chapter.  The  approximation  presented 
in  Sec.  3-3  was  first  formulated  by  Dr.  Glenn  A.  Gibson  19]. 

Table  3.1-1  is  a summary  of  the  notation  used  in  the  derivation. 
Note  that  the  Lorentz  half-width  is  given  simply  as  "a"  and  the  sub- 
script "1"  now  refers  to  the  atmospheric  level  under  consideration. 

3.2  Derivation  of  the  Transmittance  Problem 
Consider  the  radiant 

energies  1^  and  I^  as 
shown  in  Fig.  3.2-1.  Ray 
1^  is  attenuated  by  the 
molecules  in  the  gas  amount 
labeled  U^,  Ray  I.  is 
incident  at  zenith  angle 
0 wrt  (with  respect  to) 
the  z-axis  and  must  pass 


VARIABLE 


MEANING 


II. 


b 

c 

E 

g 

h 

K 

k 

M 

P 

Q 

S 

T 

U 

a 

T 

V 

P 


gas  constant 
speed  of  light 
ground-state  energy 
acceleration  due  to  gravity 
Planck's  constant 
absorption  coefficient 
Boltzmann's  constant 
gas  mixing  ratio 
atmospheric  pressure 

ratio  of  vibrational  partition  functions 

spectral  line  strength 

atmospheric  temperature 

gas  amount 

Lorentz  half-width 

transmittance 

frequency 

gas  density 


SUBSCRIPT 


IDENTIFIER 


i 

j 

1 

o 


spectral  line 
molecular  gas 
atmospheric  level 
reference  or  standard 


Table  3.1-1  NOTATION  INDEX 
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through  a gas  amount  equal  to  If  is  known,  then  the  gas 

amount  for  any  slant-path  angle  can  be  calculated  by 


U2  = (sec  0)U^ 


3.2-1 


Therefore,  Eq.  2.4-4  can  be  rewritten  to  include  an  arbitrary  slant- 


path  angle  as 


t = EXP  -sec  0 K(z)  p(z)  dz  , 
o 

Z1 


3.2-2 


where  dz  is  the  vertical  path  length,  and  0 is  measured  as  shown  in 
Fig.  3.2-1. 

It  is  normally  more  convenient  to  write  Eq.  3.2-2  in  terms  of 
temperature  and  pressure.  The  hydrostatic  approximation,  as  given 
by  Eq.  2.1-1,  can  also  be  written  as 


”^=p(z)  dz 
g 


3. 2-2a 


where  M(P)  is  the  gas  mixing  ratio.  Using  this  form,  transmittance 


is  now 


t(v,T,P)  = EXP 


- sec  0 


K(v,T,P)  M(P)  dp  . 3.2-3 


In  order  to  arrive  at  an  iterative  solution  to  Eq.  3.2-3,  it  is 
desirable  to  express  t as  a product  of  a known  transmittance  and  a 


'loss"  function.  Hence,  in  equation  form,  we  seek 


XL+1  = TL  f(v’T’P)’ 


3.2-4 
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where  f(v,T,P)  is  the  "loss"  function,  and  is  known. 


As  Illustrated  in 


Fig.  3.2-2,  we  would 
like  to  express 
transmittance  at  any 
level,  tl+^>  ^ terms 
of  the  previous  known 
value,  tl»  To  do  this 
we  first  write  t,  as 


1 = first 


1 « last 


Figure  3.2-2 


defined  in  Eq.  3.2-3,  for  two  successive  levels.  Hence,  for  the  Lth 
level 


t * EXP 
Li 


-sec  0 


K(v,T,P)  M(P)  dP  , 3 . 2-5a 


first 


and  for  the  (L+l)th  level. 


Vl  " EXP 


-sec  6 


K(v,T,P)  M(P)  dp 


first 


Noting  that  pfirst  <_  i.  pL+i  and  using  a well-known  property  of 

0 

definite  integrals  [Fuller  ] , can  t>®  expressed  as  the  sum  of  two 

integrals  as 


TL+1  " ^ ~gC  I j “M*  + JpL+1  KMdP  j , 3 . 2- 5b 

Pf irst  L ' 


where  the  terms  K.(v,T,P)  and  M(P)  have  been  abbreviated.  Now  we 


divide  Eq.  3.2-5b  by  Eq.  3.2-5a  and  solve  for  tl+^.  This  gives 


TL+1  = TL 


-sec  0 


3.2-6 


If  we  define  the  "loss"  function  as  the  term  in  brackets,  the  desired 
form  of  Eq.  3.2-4  is  achieved. 

A number  of  techniques  exist  to  solve  Eq.  3.2-6  in  its  integral 
form  for  ICarnahan^].  The  expression,  however,  can  also  be  trans- 

formed into  a differential  equation  which  may  be  easier  to  solve  numer- 


ically. 


Taking  the  natural  logarithm  of  Eq.  3.2-6  we  get 


ln(xL+l)  “ ln(TL)  + 


-sec  0 


KMdp 


Since  we  assume  t.  is  known,  ln(x  ) is  a constant  and  its  derivative 
L L 

will  be  zero.  Therefore,  differentiating  wrt  pressure  gives 


1 L+l 

tl+i  dp 


-sec  0 


From  elementary  calculus  we  know  that 


d If  f(x)  dx] 

—*5 £W' 


Hence,  the  above  equation  can  be  written  as 


"IT1  ’ tl+i  [ T5"  KM  ] ‘ 3-2-7 

Equation  3.2-7  is  a first  order,  ordinary  differential  equation 
with  initial  conditions  that  at  P1  = P^,  = t^.  Many  standard  solu- 

tions exist  for  the  equation  in  this  form  [Carnahan^]  and  it  is  common- 
ly believed  to  be  easier  to  computerize. 

To  solve  for  from  either  equation  derived  above,  it  is  necessary 

to  know  several  quantities.  They  are:  (1)  slant-path  angle  - 0;  (2) 
gas  mixing  ratio  - M(P) ; and  (3)  absorption  coefficient  - K(v,T,P). 

Problem  geometry  will  always  dictate  angle  0 and  M(P)  can  be  found 
easily  from  techniques  discussed  in  Chapter  4.  Therefore,  for  the  pur- 
pose of  this  chapter,  K(v,T,P)  is  the  only  remaining  unknown.  Section 
3.3  will  formulate  a method  of  calculating  the  absorption  coefficient 
assuming  a Lorentz  line  shape. 


3.3  Derivation  of  the  Approximation 

In  the  pressure  regions  from  1013  mb  to  100  mb,  spectral  line 
broadening  may  be  accurately  described  by  the  Lorentz  contour  [McClat- 
chey^].  Using  the  notation  of  Table  3.1-1,  K(v,T,P)  can  be  written 
for  the  jth  absorbing  gas,  at  a given  level,  as  the  summation  over  all 
contributing  lines: 


Kj(v,T,P) 


fil 

it  2 , 

“u + (v-w 


3.3-1 
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Additionally  the  line  strength  and  half-width  can  be  written  as 


3.3-2 


3.3-3 


Assuming  that  a temperature  and  pressure  profile  Is  known  snd 
the  line  parameter  data  for  ^oij  * a0ij  » s0ij  ancI  Eij  *s  available,  we 
wish  to  calculate: 

(1)  t(v,T,P)  for  a fixed  frequency  vq 
or 

(2)  x(Av,T,P)  for  a band  of  frequencies  Av. 


Regardless  of  which  transmittance  Is 
time  will  be  spent  on  evaluating  Gq. 
duced  here  will  shorten  this  time  by 
ations  that  need  computing  only  once 
tral  line. 

With  this  aim  in  mind,  we  let 


desired,  the  bulk  of  calculation 
3.3-1.  The  approximation  intro- 
separating  those  arithmetic  oper- 
for  a given  level,  gas,  or  spec- 
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where,  for  simplicity,  the  "j"  subscript  has  been  omitted  as  all 
calculations  refer  to  the  jth.  absorbing  gas.  From  these  definitions, 
Eq.  3.3-1  becomes 


K(v,T,P) 


- J noi  . v2  , 2 • 

i oi  Xi  + Yi 


3.3-4 


The  expression  —z can  be  broken  up  into  once-per-level , and 


Xl  + Yi 


once-per-line  calculations.  For  an  arbitrary  positive  constant  "a". 


2 2 


2 2 
Y*  + a 


2 2 
3 ~X1 

Y2  + a2 


3.3-5 


2 2 
a “ X1 

If  2 2 < C^en  3.3-5  can  be  expanded  by  a binomial 

Y*  + a 


series  [Selby  ] as 


h . xi  V / ‘2-xi 

+ Yi  Yi + »2  \ ?i + *2 


Defining 


fmCYi)  " (Yi  + a2)“(m+1)  and 


*„«!>  - XjCa2  - xj)  , 


Eq.  3.3-4  is  now 


KCv.T.P) 


V~  V si 

L.  L,  ^8»(xi,£»(Yi)' 

i.  m*»0 


3.3-6 


The  next  expression  to  be  separated  is  found  in  Eq.  3.3-2  as 


EXP 


i 


Let  **  — , Eg  and  be  constants  chosen  to  be  in  the  middle  of 
their  respective  ranges  of  interest.  Then 


44  t1-^) 


EXP 


E E.  - E 1 
a i a 

kT  kT  l 

o o > 


1 - W - W,  + W 
ala 


Noting  that  E*  >=  ^ "jy  ISelby^],  then 


E T 

EXP  ~ 1 - -2. 

kT  T 

o 1 


E -E. 

-£t-J  W1  " Wa 


Using  this  expression,  S (from  Eq.  3.3-2)  is 


Si  " W 


'Wo,  t 


[WIK 


3.3-7 


Finally,  using  Eqs.  3.3-6  and  3.3-7,  the  Lorentz  line  absorption 
coefficient  for  the  1th  level  and  jth  gas  is  exactly 


1 bl  (W,  - w )n  “ 

K«  " 7 % *i  "u  Z — in-5—  Z 8»  (xi) 


? I 


E . E. 


V1!  W 


3.3-8 


Equation  3.3-8  is  the  basis  for  absorption  coefficient  calculations 

I 

in  which  quantities  dependent  on  atmospheric  level,  molecular  gas,  and 
spectral  line  have  been  separated. 


'I 


3 . 4 Selection  of  Arbitrary  Constants  and  Series  Truncation 

Three  decisions  involving  constants  must  be  made  before  Eq.  3.3-8 
can  be  used  to  approximate  the  Lorentz  line  absorption  coefficient. 
They  are: 

(1)  selecting  appropriate  values  for  E and  W ; 

fl  a 

(2)  fixing  the  constant  "a"  in  Eq.  3.3-5  to  insure 
convergence  of  the  binomial  series;  and 

(3)  truncating  the  infinite  series  over  "m"  and  "n" 
so  that  a desired  degree  of  accuracy  is  achieved 
with  the  minimum  number  of  terms . 

The  line  parameter  data  used  in  this  document  [McClatchey^ ] 
was  compiled  at  a standard  temperature  (Tq)  of  296°  K.  Using  Tq  - 296° 
K and  choosing  T to  be  an  average  temperature  selected  from  the  pro- 

d 

file  used,  we  can  pick  W : 

Si 


w 


296 

T 


3.4-1 


E should  also  use  an  average  value  energy  for  a given  gas. 
aj 

For  the  user's  convenience,  average  energies  are  computed  and  furnished 
by  the  computer  program  (see  Appendix  A)  used  to  extract  data  from  the 
line  parameter  tape.  These  values  of  E are  used  in  all  computations 
throughout  this  work. 

Secondly,  the  expansion  of  Eq.  3.3-5  requires  the  constraint  that 


a2-x2 

y2+a2 


< 1 


3.4-2 
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1 


for  all  real  values  of  x and  y.  For  typical  atmospheric  profiles 
from  0 to  100  km  [McClatchey32] , the  ranges  of  temperature  and  pres- 
sure are  approximately  300°  K to  200°  K and  1.4  x 103  mb  to  5.0  x 10~2 
mb  respectively.  Consequently,  if  Tq  = 296°  K and  Pq  = 1013  mb  are 
taken  to  be  the  standard,  x will  range  from  near  zero  to  about  1.3, 
and  Y can  vary,  theoretically,  from  0 to  Using  these  bounds  for 
x and  y we  can  now  choose  an  "a"  and  analyze  the  constraint  for  regions 
of  convergence  and  divergence. 

Inequality  3.4-2  requires  two  conditions  to  be  met  depending  on 

the  sign  of/  2_x2^  When  a>x,  the  inequality  reduces  to 

\~2~2 

\y  +a 


2 2 
y >-x 


2 2 

Since  x and  y are  real,  -x  will  always  be  negative  and  less  than  y . 

Hence,  the  series  will  always  converge  for  a>x.  Next,  when  a<x  the 

Inequality  becomes 

2 2 2 

y < x - 2a  • 3.4-3 

This  restriction  will  be  most  severe  when  y = 0.  When  this  occurs, 
inequality  3.4-3  becomes 

x2  < 2a2  • 3.4-4 

Inequality  3.4-3  is  plotted  in  Fig.  3.4-1  using  only  the  I and  IV 
quadrants  (regions  of  interest)  for  two  values  of  "a".  Based  on  Eq. 
3.4-4,  convergence  can  always  be  assured  when 
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a > 


max 

/2~ 


3.4-5 


where  is  the  upper  bound  of  x.  Satisfying  Eq.  3.4-5  will  push 

the  region  of  divergence  from  the  area  of  valid  x and  y values.  This 
is  illustrated  in  Fig.  3.4-1. 

Insuring  the  convergence  of  the  binomial  series  and  selecting  W 

a 

and  E have  been  shown  to  be  straight  forward.  Truncating  the  two 
aJ 

infinite  series  in  Eq.  3.3-8,  however,  requires  a more  complicated 
approach.  Because  a finite  number  of  terms  must  be  used,  an  error  is 
introduced  into  the  approximation.  The  magnitude  of  the  error  can  be 
made  arbitrarily  small  by  increasing  the  number  of  terms,  but  this  has 
the  adverse  effect  of  increasing  calculation  time.  Care  must  be  taken 
in  keeping  the  number  of  terms  to  a minimum  without  allowing  the  result- 
ant error  to  become  larger  than  desired. 

First,  consider  the  Taylor's  series  expansion  used  in  Eq.  3.3-7 
which  is  repeated  here: 


EXP 


IvM 

(urBa ) 

00 

- Y ±- 

’(E.-V  | 

[w.-w  ' 

1 a 

n 

. 3.4-6 

l kT  1 

Zj  n! 

kT 

. 0 1 

mJ 

n-0 

0 

tm 

1 J 

Let  Z = 

J kT„  1 * 

wrwa] 

. The  average  error  after  n+1  terms  is 


given  as 


R(Z) 


Za 


n+1 


exp  [z  3, 

(n+1) 1 max 
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is  the  maximum  value 


where  Z is  the  average  value  of  Z and  Z 

a max 

I Carnahan^ ] . Using  typical  values  from  line  parameters  in  the  4.3pm 
band  for  channels  11--15  we  have, 


'Wave  ; 2201-°°“ 


Wave  4 °'199 


'Wmax  ' 2381 ' 602 


WW  ■ °-222 


(see  McClatchey  [17]  and  Table  3.4-1).  If  the  units  of  (E^-E^)  are 

-4 

converted  properly  and  if  we  choose  R(Z)  = 5.0  x 10  to  insure  three 
decimal  place  accuracy,  we  establish  the  error  bound  as 


5.0  x 10~4  _>  fl*0]  EXP  £ 2.500  J 


which  is  satisfied  when  n = 12.  A further  reduction  in  the  number  of 
terms  needed  can  be  realized  through  Chebyshev  economization  of  EXP[Z], 
A reduction  to  nine  terms  was  obtained  in  the  4.3pm  band. 

Lastly,  the  binomial  series  of  Eq.  3.3-5  must  also  be  truncated. 
Recall  that 


2 2 
Y + X 


v?  iM 


3.4-7 


.2  ,2  \“ 

Let  V = --A-—  \ j . Then 

, Y + a 


w 

-j_.y 

i - v Z_i 


m+l 

vrn  = i + v + v2+  ’ ’ * + V + f— y , 


3.4-8 


1148.477 

710.910 


1100.600 

1546.528 
1469.410 
522.473  I 224.378 


661.560  1219.815 

27.359  2381.602 

37.710  1442.873 

1346.596  401.388 


446.710  I 1348.960 
2293.701 

653.221  | 1164.6845 


2549.953 


0.000 


542.910  I 1385.995 


294.057 


2334.756 

562.775 


*Units  of  E are  in  cm 


Table  3.4-1  GROUND  STATE  ENERGY  DISTRIBUTION  IN  THE  4.3ym  BAND 
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and  the  exact  error  is 


V 


R = 


yBri-l- 
1 - V 


3.4-9 


From  Eq.  3.4-9,  the  absolute  fractional  error  can  be  formed  : 


R 


f 


vmn 

1 - V 

vm+1 

/ a2  - X2 

1 

1 - V 

U+a2j 

Gibson  [10]  has  determined  that  for  a fractional  error  of  one  part 
in  ten  thousand,  the  optimum  value  for  the  constant  "a"  will  not  com- 
pletely satisfy  the  convergence  criteria  of  inequality  3.4-5.  Diver- 
gence of  the  binomial  series  takes  place  whenever  | (v-vo) | <0.4 
wavenumber.  Since  this  occurs  often,  the  series  of  Eq.  3.4-8  has  to 
be  modified  to  reduce  the  number  of  instances  in  which  nonconvergence 
develops. 

In  addition  to  the  value  defined  as  V,  we  add  Z: 


Then  from  Eq.  3.4-7 


X2  + Y2  1 ~ V 


] 


2 =z  i + v + v2+...+vm  + 


ft] 


,nrt-l 


If  the  last  term  in  series,  ^ - , is  approximated  by 
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cv**1  = 


. v"*1 


1 - V 


where,  for  an  arbitrary  constant  "d" , 


d2  (Y2  + a2) 


(X2  + d2)  Y2 


For  M+l  terms, 


2 2 
X + Y 


= Z 


H-l 

I * 


m=0 


ZV^C 


M-l 


m»0 


X(a2  - X2)m 
(Y2  + a2) 111+1 


2 2 
d X(a 


x2)M 


2 2 
(X  + dZ) 


Y2(Y2  + a2)M 


3.4 


Equation  3.4-10  has  the  effect  of  adding  the  correction  factor 
Z^C  to  the  approximation  of  Eq.  3.4-8.  The  error  term  can  now  be 
reduced  by  selection  of  constants  "a"  and  "d".  Unpublished  analysis 
[Gibson1^]  suggests  that,  using  7 terms  and  the  same  fractional  error 
of  .0001,  divergence  of  Eq.  3.4-10  can  be  avoided  unless  | (v  - v ) | < 
0.0525  wavenumber.  Fortunately,  this  will  happen  in  less  than  0.5%  of 
the  spectral  lines  considered.  For  these  few  cases  it  is  not  incon- 
venient to  calculate  the  contribution  to  the  absorption  coefficient 
using  the  standard  formula  of  Eq.  3.3-1. 

It  will  be  shown  in  Chapter  6 that  calculation  of  transmittance 


-10 
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using  Eq.  3.3-8  and  the  modification  explained  above  is  accurate  to 
at  least  three  decimal  places.  Additionally,  it  will  be  seen  that 
because  of  the  separation  technique  used  in  the  approximation  and 
the  small  number  of  terms  required  for  series  convergence,  a signifi- 
cant amount  of  computational  time  is  saved. 
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CHAPTER  IV 
i" 

TOPICS  RELATED  TO  ATMOSPHERIC  ABSORPTION 

4.1  Preface 

This  chapter  is  concerned  with  developing  the  tools  necessary 

to  apply  the  theory  presented  in  Chapters  2 and  3.  Before  Eq.  3.3-8 

can  be  used  to  calculate  the  Lorentz  line  absorption  coefficient, 

spectral  line  parameters,  v , a , S and  E"  must  be  known.  Addi- 

o o o 

tionally,  the  various  formulae  will  yield  meaningful  results  only 
if  proper  units  and  conversion  factors  are  employed.  Because  a 
working  knowledge  of  units  is  essential  to  the  atmospheric  investi- 
gator and  the  amount  of  line  parameter  data  one  must  handle  can  be 
overwhelming,  this  chapter  will  address  some  of  the  practical 
problems  often  encountered  while  working  in  atmospheric  sciences. 

4 . 2 Line  Parameters  - AFCRL's  Data  Tape 

An  excellent  source  of  spectroscopic  data  is  the  Air  Force 
Cambridge  Research  Laboratories'  (AFCRL)  Atmospheric  Absorption  Line 
Parameters  Compilation  [McClatchey ^ ] . The  tape  includes  data  on 
108,000  spectral  lines  ranging  in  wavelength  from  . 828pm  to  .0325m. 
Seven  molecules  are  included  in  the  compilation  which  are  found 
naturally  in  the  earth's  atmosphere  and  considered  to  be  the  absorbers 
in  the  infrared  and  visible  spectrums.  They  are:  (1)  water  vapor  - 
^0;  (2)  carbon  dioxide  - CO^;  (3)  ozone  - 0^;  (4)  nitrous  oxide  - 
N2O;  (5)  carbon  monoxide  - CO;  (6)  methane  - CH^ ; and  (7)  oxygen  - 


43 
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To  enable  the  user  to  implement  the  equations  of  Chapters  2 

and  3,  four  essential  parameters  are  included  for  each  data  entry 

along  with  the  molecule's  identification  information.  A spectral 

line  dossier  includes:  the  resonant  frequency,  vq  (cm  ^) ; the  inten- 

-1  -2 

sity  per  absorbing  molecule,  So  (cm  /molecule-cm  );  the  Lorentz 
half-width,  aQ  (cm  ^);  the  rest-state  energy,  E"  (cm  ^) ; the  rota- 
tion and  vibration  identification  numbers;  date  of  information;  the 
isotope  identification  code;  and  lastly,  the  molecule  identification 
number.  The  lines  are  ordered  on  the  tape  in  ascending  frequency 
(measured  in  wavenumbers) . 

Lines  were  included  in  the  compilation  based  on  an  intensity 
criteria.  The  "Criterion  Intensity  Minimum"  for  each  molecule  is 
shown  in  Table  4.2-1.  The  cutoff  values  were  violated  in  two 


Molecule  Type 

I.  D.  No. 

S Minimum>' 
o 

H2° 

1 

-27 

3.0  x 10 

co2 

2 

2.2  x 10~26 

°3 

3 

•3  c in~24 

3.5  x 10 

N20 

4 

-23 

3.0  x 10 

CO 

5 

-23 

8.3  x 10 

CH. 

4 

6 

3.3  x 10~24 

°2 

7 

7 in-30 

3.7  x 10 

-1  -2 
*Units  are  cm  /molecule-cm 

Table  4.2-1  CRITERION  INTENSITY  MINIMUM 
[Adapted  from  McClatchey^] 


instances.  In  regions  of  very  strong  absorption,  very  weak  lines 
above  the  established  minimum  were  omitted.  Also  for  cases  in  which 
weak  lines  are  so  closely  spaced  that  their  combined  effect  may 
contribute  significantly  to  absorption,  the  minimum  criteria  was 
lowered . 

Table  4.2-2  is  a listing  of  the  coding  used  to  identify  each 
molecule  type  and  isotope.  Molecules  are  identified  by  a simple 
numbering  system  from  1 to  7.  Isotope  codes  are  based  on  the  mole- 
cule's atomic  numbers.  For  example,  the  four  isotopes  of  carbon 
monoxide  can  be  written  as:  CO;  CO;  CO;  CO.  If  the 

units  digit  of  the  atomic  numbers  are  used,  the  respective  isotope 
codes  for  carbon  monoxide  are:  26;  36;  28;  and  27. 

The  spectral  data  compiled  by  AFCRL  is  on  a 7-track,  800  BPI, 
1/2  in.  wide  magnetic  tape.  The  format  used  for  each  line  entry  is 
standardized  to  an  80  column  card  and  is  given  below.  Data  records 


Variable 

V 

0 

s 

o 

a 

o 

E" 

Rot.  & Vib. 

I.  D.  No. 

Date 

Isotope 

Molecule 

Card  Column 

1-10 

11-20 

21-25 

26-35 

36-70 

71-73 

74-77 

78-80 

Format  Code 

F10.3 

E10.3 

F5.3 

F10.3 

5 A6,  A5 

13 

14 

Ie 

Figure  4.2-1  WORD  FORMAT  CHART 

are  composed  of  a number  of  card-imaged  line  parameter  sets.  Each 
record  begins  with  a number  (format  110)  indicating  how  many  card- 
imaged  sets  are  to  follow  and  is  separated  by  a string  of  alfanumeric 
characters  and  symbols  which  aid  in  record  identification.  A full 


■■■■ 
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Table  4.2-2  MOLECULE  AND  ISOTOPE  CODES 
[Adapted  from  McClatchey^] 
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record  contains  40  card  images  in  the  above  format. 


End  of  file  (EOF) 


markers  are  placed  on  the  tape  at  frequency  points:  500,  1000,  2000, 
5000,  7500,  and  10000  cm  A double  EOF  marker  is  at  the  end  of  the 
tape. 

4.3  Extracting  Data  from  AFCRL's  Tape 

In  the  report  referenced  above  [McClatchey ^ ] , a tape-read 
program  is  listed  which  can  be  used  with  a Control  Data  Corporation 
660  computer.  Since  there  is  a good  chance  that  some  users  may  not 
have  CDC  computer  equipment,  a data  retrieval  routine  is  listed  in 
Appendix  A for  use  with  an  IBM-360/65. 

The  routine  given  in  the  Appendix  was  written  to  perform 
several  tasks: 

(1)  access  AFCRL's  tape  and  extract  spectral  line  data 
between  arbitrary  frequency  limits  (channel  limits) ; 

(2)  count  the  number  of  lines  contained  within  the  channel 
for  each  molecule  type; 

(3)  segregate  line  parameters  by  molecule  type  into  separate 
arrays ; 

(4)  build  a single  data  file  for  arbitrary  frequency  limits 
with  data  arranged  in  ascending  frequency  and  segregated 
by  molecule  type;  and 

(5)  compute  the  mid-point  value  of  energy  and  the  average 
intensity  for  each  molecule  type. 

The  program  is  divided  into  two  parts.  The  main  routine  Initializes 
variables,  selects  channel  limits  and  outputs  all  data.  The  subroutine. 
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CHREAD,  reads  the  tape  between  channel  limits,  counts  the  number  of 
lines  for  each  molecule  type,  segregates  the  data,  writes  channel 
data  on  a permanent  file  and  computes  average  values  of  energy  and 
intensity  described  in  (5)  above.  The  user  can  easily  adapt  all  or 
part  of  the  routine  for  a particular  purpose.  Flow  charts  for  the 
main  and  subroutine  programs  can  be  found  in  Fig.  4.3-1.  Table  4.3-1 
is  a listing  of  the  main  variables  used  in  the  programs. 

4.4  Units  of  Atmospheric  Physics 

This  section  is  presented  to  the  novice  atmospheric  scientist 
as  an  aid  in  understanding  the  units  used  in  infrared  absorption 
calculations.  For  the  sake  of  continuity,  the  subject  of  units  will 
be  discussed  as  it  applies  to  the  basic  equation  of  transmittance. 
First,  however,  a clarifying  note  is  in  order  to  explain  the  term 
"frequency"  which  is  used  rather  loosely  in  atmospheric  physics. 

In  almost  all  cases,  frequency  will  be  synonymous  with  the  term 
wavenumber  (cm  ^) . The  wavelength,  X,  and  wavenumber,  v,  are  related 

by  4 

, -1,  10 

Vlcm  ] = TUm]  ’ 


where  X is  given  in  micrometers.  For  example. 


Variable 


Explanation  of  Variable  Name 


AN2Q 


Data  array  for  nitrous  oxide 


BWDTH 


Array  containing  channel  bandwidths 


CH4 


Data  array  for  methane 


CNTR 


Array  containing  channel  center  frequencies 


CO 


Data  array  for  carbon  monoxide 


C02 


Data  array  for  carbon  dioxide 


EAVE 


Array  containing  Eafor  each  molecule  type 


FHI 


Channel  upper  frequency 


FLO 


Channel  lower  frequency 


H2O 


Data  array  for  water  vapor 


ICNT 


Array  containing  number  of  lines  per  molecule 


IFILE 

KMOL 

NOFIL 

02 

03 

SAVE 


found  within  channel  limits 

Array  used  to  reference  a data  file 

Molecule  identification  number 

Data  file  reference  number 

Data  array  for  oxygen 

Data  array  for  ozone 

Array  containing  average  intensities 


Table  4.3-1 


VARIABLE  NAMES  FOR  DATA  RETRIEVAL  ROUTINE 
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SEGREGATE  CARD-IMAGE  DATA 


SAVE  MAX  AND  MIN  VALUES  OF 
ENERGY  AND  SUM  INTENSITY 


GO  TO  FILE  READ  BLOCK 


OUTPUT  LINE  COUNTS 
FOR  CHANNEL 


Figure  4.3-l(cont.)  CHREAD  DATA  ROUTINE  FLOW  DIAGRAM 


51 


Figure  4.3-l(cont.)  CHREAD  DATA  ROUTINE  FLOW  DIAGRAM 
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The  dimensions  of  the  line  parameters  (a  , S , E",  v ),  temperature, 

o o o 

and  pressure  are  normally  established  by  the  data  available  for  use. 
Typical  dimensions  for  these  values  are  shown  below. 
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Variable 

Dimension 

T 

®K  (degrees  Kelvin) 

P 

mb  (millibars) 

a 

o 

cm  '*'/atm  (wavenumber/atmosphere) 

-1 

V 

o 

cm  (wavenumber) 

S 

o 

cm  ^ /molecule  - cm  ^ 

-1 

E" 

cm  (wavenumber) 

Pressure  is  also  commonly  given  in  atmospheres  [atm] , and  is 
related  to  millibars  [mb]  by 

l[mb]  = 9.869233  x lCf4  [atm].  4.4-5 

Since  the  standard  pressure,  Pq,  used  throughout  AFCRL's  data  com- 
pilation is  one  atmosphere  (P  = latm  - 1013  mb) , the  numerical 

value  of  the  half-width  is  actually  a /P  dimensionally.  Therefore, 

o o 

the  dimensions  of  aQ  are  really  [cm  ^]  and  will  be  used  as  such. 
Lastly,  the  partition  function  ratio,  Q,  and  Boltzmann's  constant, 
k,  which  appear  in  Eq.  4.4-4,  have  dimensions  of 


and 

k = Joule/°K, 

where  « means  "dimensionally  equal  to." 

With  these  facts  in  mind,  Eqs.  4.4-3  and  4.4-4  can  be  written 
dimensionally  as 
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d r ' 
a ■ Icm 


‘{iH 


1/2 


d 

3 


4.4-5 


and 

e d 

-1 

cm 

Til 

°K 

b 

EXP 

cm-1  °K 

b =* 

molecule-cm  2 

rJ 

°K 

(Joule/°K) (°K2) 

P — 

— — 

-1 

-1 

cm 

EXP 

cm 

molecule-cm  _ 

Joule 

4 .4-6a 


The  dimensions  of  the  exponential  in  Eq.  4.4-1  must  be  made 
equal  to  unity.  This  can  be  accomplished  by  converting  ground-state 
energy  from  wavenumbers  to  Joules.  First,  note  that  Planck's 
constant,  h,  and  the  speed  of  light  in  a vacuum,  c,  are  given  as 
[Abramowitz*] 

h = 6.626196  x 10  [Joule-sec] 

and 

c = 2.997925  x 1010[cm/sec] . 


If  the  ground-state  energy,  E",  is  multiplied  by  the  factor  he,  the 
dimensions  of  the  expression  become 
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which  is  the  desired  result.  Therefore,  energies  which  are  given  in 
wavenumbers  can  always  be  converted  to  Joules  through  multiplication 
with  "he." 

With  the  exponential  expression  dimensionally  equal  to  unity , 

Eq.  4.4-6  is  finally  given  as 

4.4-6b 

Now,  for  a and  v measured  in  wavenumbers,  the  Lorentz  line  shape 
given  by  Eq.  2.3-7  is 


■ [cm]. 


4.4-7 


Hence,  the  absorption  coefficient  is 
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molecule-cm 


-2 


4.4-8 


Equation  4.4-8  gives  K in  terms  of  units  using  only  AFCRL's 

data  tape  and  modifying  E"  as  described  above.  Converting  these 

units  to  inverse  atmospheric  centimeters  can  be  accomplished  using 

18 

elementary  chemistry  [McClatchey  ] . 

An  [atm-cm]  of  gas  at  STP  (standard  temperature  and  pressure) 

2 

is  equivalent  to  a length  of  1 cm  of  gas  at  STP  per  cm  . From 

Avogadros'  hypothesis,  the  molecular  weight  (Mt) , in  grams,  of  any 

-3  3 

gas  at  STP  occupies  a volume  of  22.4  x 10  cm  (22.4  liters). 
Therefore, 

Mt  2 

1 [atm-cm]  = r [gm/cm  ], 

2.24  x 10 


at  STP.  From  the  fact  that  one  mole  of  gas  contains  the  same  number 
of  molecules  (Avogadro's  constant)  and  one  cubic  centimeter  would 


contain 


2.24  x 10H 


moles , it  can  be  deduced  that 


1 [atm-cm] 


STP 


.6884683  x 101 ^ [molecules/cm^] . 


4.4-9 


Therefore,  K in  Eq.  4.4-8,  can  be  given  in  inverse  atmospheric 
centimeters  using  Eq.  4.4-9: 

K » C[atm-cm]  \ 4.4-10 

where  C is  the  conversion  factor  defined  above. 
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In  summary,  it  is  relatively  simple  to  express  K in  inverse 
atmospheric  centimeters.  First  E"  must  be  converted  to  Joules 
before  the  line  strength  equation  can  be  used.  Secondly,  K must  be 
multiplied  by  the  conversion  factor  defined  in  Eq.  4.4-9.  Other- 
wise, for  a given  atmospheric  profile,  the  AFCRL  data  tape  can  be 
used  directly  with  Eqs.  4.4-2,  4.4-3,  and  4.4-4  to  calculate  K. 

The  gas  amount,  U,  must  now  be  expressed  in  terms  of  atmospheric 
centimeters.  From  Eq.  2.4-1,  the  differential  gas  amount  can  be 
expressed  as 

p (z) 

du  5 — ° dz  [atm-cm],  4.4-12 

PSTP 


where  p is  the  density  of  the  gas  at  STP.  Equation  4.4-12  is  not 

bit 

commonly  used,  however,  since  parameters  are  normally  given  in  terms 
of  pressure  profiles  and  gas  mixing  ratios . 

If  the  mixing  ratio  for  a gas  is  defined  in  terms  of  its  density 
to  that  of  dry  air,  we  have 

M = . 

pair 


air 


Using  this  definition  and  multiplying  Eq.  4.4-12  by  , we  have 


air 


M p 


dU 


air 


STP 


4.4-13a 


From  the  hydrostatic  equation,  the  altitude,  z,  can  be  expressed  in 
terms  of  pressure  so  that  the  gas  amount  is  now 
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4.4-13b 


dU  - 


103-M-dP 
g PSTp 


P is  given  in  JmbJ;  g is  the  acceleration  due  to  gravity  in 
2 

[cm/sec  ];  and  M is  dimensionless. 

Equation  4.4-13b  is  used  for  gases  which  have  a variable  mixing 
ratio  (i.e.  1^0  and  0^) . Using  the  definition  of  M,  if  the  densi- 
ties of  the  gas  and  of  dry  air  are  known,  a mixing  ratio  profile  can 
be  computed  easily. 

For  the  constant  mixing  ratio  gases  (see  Sec.  2.5),  the  gas 
amount  is  given  by 


dU 


V”,,-10'3  jp 

Mta’8 *'pSTP 


4.4-14 


where  Mpp^  is  the  "parts  per  million"  mixing  ratio,  and  M and  Mta 
are  the  molecular  weights  of  the  gas  and  air  respectively.  Since 
these  quantities  are  easily  found,  Eq.  4.4-14  can  be  written  in 
terms  of  a constant  times  a pressure  differential  for  a given  gas: 


dU  * A dP, 
g 


where 


“ppm’11, 


10 


-3 


8 "“ta "PSTP 


Table  4.4-1  is  a listing  of  molecular  weights,  constant  parts 

per  million  mixing  ratios  and  gas  densities  at  STP.  Also  listed  are 

computed  values  of  A where  appropriate. 
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Mixing 

Ratio-PFM* 

Gas  Density 

At  STP** 

A Constant 

8 

of  Eq.  4.4-15 

h2o 

18 

Variable 

1.00016 

N/A 

co2 

44 

330 

1.9768441  x 10-3 

2.58538380  x 10-1 

°3 

48 

Variable 

2.1436282  x 10_3 

N/A 

n2o 

44 

0.28 

1.978137  x 10~3 

2.19222523  x lO-4 

CO 

28 

0.075 

1.25101  x 10-3 

5.90866466  x 10~5 

ch4 

16 

1.6 

7.1678376  x lO-3 

1.25713708  x 10~3 

°2 

32 

2.095  x 10"5 

1.4286545  x 10_3 

1.65172454  x 10+2 

AIR 

28.97 

1 x 106 

N/A 

N/A 

18 

♦[Adapted  From  McClatchey  ] 

20 

♦♦[Adapted  From  Nelson  ] Units  In  gm/ 

3 

cm 

Table  4.4-1  GAS  AMOUNT  TABLE 
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l 


The  discussion  of  units  can  he  concluded  with  a few  remarks 


concerning  the  approximation  to  K given  in  Chapter  3.  Like  the 


standard,  Eq.  3.3-8  will  be  in  units  of 
be  converted  to  Iatm-cm]  ^ by  Eq.  4.4-10. 
Ea  must  also  be  expressed  in  iJoules], 


_1 

. -2 
lle-cm  _ 


and  must 


molecule- 
Additionally  the  constant 


4.5  The  Nimbus-6  Sounder  And  The  4.3pm  Band 

When  the  Lorentz  line  absorption  coefficient  was  developed  in 

Chapter  2,  it  was  stated  that  the  line  broadening  mechanism  was  based 

on  molecular  collisions  or  pressure.  This  suggests  that  the  best 

accuracy  in  transmittance  calculations  using  the  Lorentz  line  shape 

will  be  attained  when  atmospheric  pressures  are  relatively  high. 

Commonly,  Lorentz  line  broadening  is  used  in  pressure  regions  from 

1013  mb  (sea  level)  to  100  + mb  (17-20  KM)  [McClatchey^ ) . For  lower 

pressures  the  doppler  effect  must  also  be  taken  into  account. 

Ultimately,  the  Lorentz  line  transmittance  calculations  are  to 

be  used  with  satellite  sounding  data  for  temperature  profile  retrievals 

iKaplan^].  One  such  High  Resolution  Infrared  Radiation  Sounder 

25 

(HIRS)  is  the  Nimbus-6  experiment  package  [Sissila  ].  The  Nimbus-6 
HIRS  is  a third  generation  infrared  radiation  sounding  experiment 
used  to  monitor  the  earth's  atmosphere  from  an  orbiting  satellite. 

The  radiation  measurements  in  4.3ym  and  15ym  bands  permit  high  reso- 
lution vertical  temperature  profiles  to  be  calculated  from  the  lower 
troposphere  to  the  stratopause. 

To  gather  these  soundings,  the  Nimbus-6  has  17  frequency 
channels  from  . 69ym  to  15ym  which  are  listed  in  Table  4.5-1. 
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Channels  1-7  are  used  for  temperature  soundings  of  the  relatively 
cold  regions  of  the  atmosphere  and  for  cloud  detection.  Channels  8 
and  16  are  used  to  measure  surface  radiances.  Water  vapor  concen- 
trations are  monitored  through  channels  9 and  10.  Similar  to  chan- 
nels 1-7,  temperature  soundings  and  cloud  detection  of  the  warmer 
regions  of  the  atmosphere  are  accomplished  with  channels  11-15. 
Channel  17  is  used  primarily  for  cloud  detection  during  daylight 
hours . 

In  selecting  a frequency  region  to  test  the  approximation  to 
the  Lorentz  line  absorption  coefficient,  the  Nimbus-6  specifications 
were  studied.  Channels  11-15  were  chosen  because  the  4.3pm  band  is 
primarily  used  for  temperature  soundings  for  regions  in  which  Lorentz 
line  broadening  is  valid.  This  is  illustrated  in  Fig.  4.5-1  which 
depicts  the  weighting  function  curves  for  channels  11-15.  Except 
for  channel  15,  all  curves  are  confined  to  high  pressure  regions. 

Also  shown  is  the  corresponding  channel  filter  function  for  the 
Nimbus-6.  The  data  contained  on  AFCRL's  tape  for  the  4.3pm  band  is 
summarized  in  Table  4.5-2.  The  information  contained  in  the  table 
was  obtained  using  the  data  retrieval  program  given  in  Appendix  A. 


08  k N20 


0.0  L_ 
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Figure  4.5-1  HIRS  WEIGHTING  FUNCTIONS 
AND  FILTER  RESPONSES 
25 

[Adapted  from  Sissila  ] 


CHANNEL 
NUMBER — O 

MOLECULE  | 
TYPE  V 

11 

12 

13 

14 

15 

h2o 

34 

53 

52 

84 

42 

co2 

195 

208 

565 

1141 

432 

°3 

0 

0 

0 

0 

0 

n2o 

748 

790 

296 

101 

30 

CO 

16 

7 

8 

1 

0 

CH. 

4 

0 

0 

0 

0 

0 

°2  ' 

0 

0 

0 

0 

0 

Totals 

993 

1058 

921 

1597 

504 

Table  4.5-2  NUMBER  OF  LINES  PER  CHANNEL  BY  MOLECULE  TYPE 
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CHAPTER  y 


TESTING  THE  APPROXIMATION 


5.1  Preface 

The  standard  formula  for  evaluating  the  absorption  coeff^ient 
due  to  Lorentz  line  broadening  is 


K = 


If 

i**l 


2 ^ , ,2 
o + (v-v  ) 


5.1-1 


where,  for  a given  gas,  there  are  "I"  spectral  lines  contributing  to 
the  absorption.  The  formulas  for  line  intensity,  S,  and  half-width, 
o,  are  rewritten  here  for  convenience: 


The  approximation  to  Eq.  5.1-1  was  developed  in  chapter  3 (Eq.  3.3-8). 
If  the  infinite  series  are  cut  off  after  (N+l)  and  (M+l)  terms,  then 
the  approximation  can  be  written  as 


K 


5^ 


W n 


(w-w  ) 

a_ 

nl 


I 


8m(X)* 


5.1-2 
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2 


s 

o 


E -E" 
_a 

kT 


f Or)  , 

m 


where  "I"  is  the  same  as  in  Eq.  5.1-1  and  variables  Q,  W,  n,  3.  gm» 

f , X and  Y are  defined  in  Chapter  3. 
m 

This  chapter  will  discuss  a method  of  computerizing  Eq.  5.1-2 
for  monochromatic  radiation.  Additionally,  several  methods  of  mea*- 
suring  computational  times  will  be  presented  along  with  accuracy  and 
error  prediction  techniques. 


5.2  Data  File  Description 

In  Chapter  4 the  data  retrieval  program  listed  in  Appendix  A was 
discussed.  Not  only  does  the  program  read  AFCRL's  data  tape  between 
arbitrary  frequency  limits,  but  also  it  creates  a permanent  data  file 
containing  line  parameters  \>q,  aQ,  Sq  and  E"  segregated  by  gas  type. 
By  knowing  the  number  of  data  records  per  gas  that  are  contained  on  a 
channel's  file,  the  line  parameters  can  easily  be  retrieved  for  use 
in  calculating  either  Eq.  5.1-1  or  Eq.  5.1-2.  The  "line  density" 
counts  for  channels  11-15  are  listed  in  Table  4.5-2  and  were  provided 
by  the  subroutine  "CHREAD"  (Appendix  A) . 

For  example,  channel  11  has  a center  wavenumber  of  2190  cm  ^ and 
a bandwidth  of  20.6  cm  ^ with  "line  densities"  of:  (1)  1^0  with  34 
lines;  (2)  CO^  with  195  lines;  (3)  ^0  with  748  lines;  (4)  CO  with 
16  lines;  (5)  0^,  CH^,  and  02  with  0 lines.  Hence,  it  is  easily 

recognized  that  record  number  34  is  the  last  line  entry  for  ^0  and 
record  number  35  is  the  first  entry  for  COj. 


67 


Data  files  are  created  in  this  manner  for  the  following  program- 
ming considerations: 

(1)  Only  one  data  file  is  needed  for  a given  channel;  and 

(2)  gases  are  segregated  by  type  which  eliminates  the  need 
to  store  a molecule  identification  code,  thus  keeping 
required  storage  space  to  a minimum. 


5.3  Data  Reorganization 


Chapter  3's  development  of  Eq.  5.1-2  stressed  that  many  calcula- 
tions needed  to  compute  K could  be  done  once  and  would  never  require 
repeating  regardless  of  the  atmospheric  level  or  the  molecular  gas 


under  consideration.  In  Eq.  5.1-2,  these  quantities  appear  after  the 

summation  over  spectral  lines  I I I . For  monochromatic  radiation, 

\i=l| 

the  expression 


f W 


can  be  evaluated  once  and  stored  in  an  (N+l)  x (M+l)  array  for  each 
gas.  Therefore,  for  all  seven  gases,  only  seven  such  arrays  would 
contain  all  the  spectral  line  data  to  compute  K for  any  number  of 
temperature  and  pressure  profiles. 

Reorganizing  the  data  in  the  form  discussed  above  enables  us  to 
"take  care  of"  a problem  mentioned  in  Chapter  3 also.  In  section  3.4 
it  was  stated  that  the  binomial  series  would  diverge  for  cases  in 
which 

|(v— v )|  < 0.0525 
1 o 
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vayenumbers.  Experience  with  AJFCRL's  data  indicates  that  this  takes 
place  in  about  0.5%  of  the  spectral  lines  for  a typical  channel.  By 
identifying  these  lines  during  the  data  reorganization  portion  of  the 
program,  their  contribution  to  K can  be  computed  in  the  standard 
fashion  and  added  in  at  a later  point  without  a significant  addition 
to  the  calculation  time. 

The  computer  routine  used  to  calculate  and  compare  Eqs.  5.1-1 
and  5.1-2  is  listed  in  Appendix  B.  The  data  reorganization  portion 
discussed  in  this  section  is  done  in  the  main  routine.  It  may  be 
prudent  to  note  here  that  the  reorganization  of  data  could  just  as 
well  be  included  in  the  tape  read  program  of  Appendix  A.  This  may 
be  a more  desirable  method  if  a complete  transmittance  program  is 
being  developed. 

5.4  Evaluation  of  K 

The  evaluation  of  K by  Eq.  5.1-1,  kSTq>  and  by  Eq.  5.1-2,  K^0D, 
can  be  done  using  the  routine  listed  in  Appendix  B.  Every  effort  to 
streamline  computation  times  in  both  methods  has  been  taken.  Some 
general  guidelines  which  can  be  used  to  accomplish  this  are: 

(1)  never  divide  by  a constant  repeatedly  if  its  inverse  can 
be  stored  and  multiplication  can  bo  performed  Instead; 

(2)  only  read  line  parameters  once  from  the  file  and  store 
them  in  core;  and 

(3)  compute  only  once  and  save  quantities  dependent  on 
atmospheric  level  or  molecular  gas. 

The  routine  in  Appendix  B is  initialized  for  the  4.3pm  band 
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channel  11.  Program  comments  specify  which  cards  must  be  changed 
to  initialize  for  a different  spectral  band.  A Plow  Diagram  for 
the  routine  is  listed  in  Fig.  5.4-1.  The  results  obtained  from 
this  program  are  discussed  in  Chapter  6. 

To  conclude  this  section,  a brief  remark  is  in  order  concerning 
the  vibrational  partition  function  ratio  Q.  In  AFCRL's  data  comi- 
lation  report,  McClatchey  I 17]  gives  the  ratio  as 

Q _ Qyo(  T) 
q Q^T)  ’ 

where  is  the  partition  function  evaluated  at  T * Tq  * 296°  K. 

Data  for  Q was  taken  from  AFCRL's  report  which  lists  values  for 
0^(1)  for  each  molecule  at  seven  discrete  temperatures  ranging  from 
175°K  to  325°K.  The  subroutine  "PART"  in  Appendix  B computes  Q 
based  on  a 2nd  degree  least-squares  polynomial.  Since  Q is  contained 
in  both  Eqs.  5.1-1  and  5.1-2,  the  relative  error  induced  by  this 
method  of  computing  the  partition  function  ratio  is  zero. 

5.5  Execution  Time  Analysis 

One  of  the  primary  reasons  for  developing  the  approximation  was 
the  hope  that  the  calculation  time  of  K could  be  significantly  reduced. 
There  are  several  methods  of  estimating  in  a quantitative  measure 
just  how  significant  of  a reduction  can  be  realized.  However,  there 
is  no  absolute  measuring  device  widely  available  for  time  analysis. 
Host  timing  "clocks"  cannot  distinguish  the  times  during  program 
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MAIN  STOPI 


Figure  5.4-1  FLOW  DIAGRAM  FOR  THE  EVALUATION  OF  K- 
MAIN  ROUTINE 
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FUNCTION  SUBPROGRAM 
TEMPK  START 


execution  in  which  the  routine  is  put  in  a "wait"  state  for  some 
reason  or  another.  Hence,  the  "execution"  times  for  two  identical 
program  runs  will  very  seldom  be  the  same,  much  less  be  an  accurate 
measure  of  the  true  execution  time. 

Even  assuming  that  the  routine  to  be  analyzed  is  the  only  one 
being  compiled  in  the  computer,  it  is  highly  doubtful  that  two 
consecutive  runs  would  now  produce  the  same  execution  times.  This 
is  due  to  the  fact  that  machine  instruction  times  can  vary  because 
the  physical  location  of  the  data  being  used  does  not  necessarily 
remain  the  same  from  run  to  run.  Therefore,  it  should  be  understood 
that  any  time  comparison  figures  presented  here  are  not  absolute  and 
are  intended  to  measure  only  in  a relative  fashion. 

One  method  that  is  normally  used  to  measure  total  computation 
time  is  to  use  the  GO-STEP  time  which  is  used  in  figuring  the  number 
of  CPU  seconds  used  by  the  routine  during  execution.  This  time  is 
taken  to  be  fairly  accurate,  especially  since  it  is  directly  related 
to  money.  Most  computer  listings  show  the  GO-STEP  execution  time  in 
the  documentation  section.  For  the  University  of  Texas  at  El  Paso's 
IBM  360/65  system,  the  time  can  be  found  by  requesting  a "MSGLEVEL  ■ 
(1,1)"  resource  and  using  the  "//GO  EXECUTION"  time  which  will  be 
listed  separately. 

The  second  method  for  timing  a program  is  to  use  "clock"  routines 
which  are  available  as  internal  software  to  the  system.  The  "clock" 
is  similar  to  the  one  which  measures  the  GO-STEP  time,  but  it  has  the 
versatility  of  being  initialized  and  read  during  any  portion  of  the 
program.  Therefore,  instead  of  measuring  total  run  times  only,  the 


programmer  can  "clock"  indiyldual  sections  of  his  routine.  Start  and 
stop  commands  can  be  obtained  through  the  computer  facility  haying 
these  internal  subroutines. 


In  case  such  a routine  is  not  available,  one  is  provided  in 
Fig.  5.5-1.  This  routine  is  written  in  Assembler  and  is  linked  to 
the  Fortran  program  by  the  Job  Control  Language  (JCL)  provided.  Out- 
put is  given  in  microseconds  (psec) . In  step  5 of  Fig.  5.5-1  an  ex- 
ample of  using  the  clock  to  time  the  subroutine  "ABSTD"  is  given.  The 
clock  is  initialized  by  calling  "STIMER"  and  when  the  program  returns 
from  the  subroutine  the  clock  is  stopped  and  read  by  using  the  "TTIMER" 
call  statement.  The  output  will  be  the  elapsed  time  in  CPU  psec  plus 
any  time  injected  while  the  routine  was  in  a "wait"  state. 

Lastly,  a different  method  which  does  not  use  a "clock"  to  com- 
pare computation  times,  is  presented  here.  A relative  time  analysis 
can  be  done  by  counting  every  arithmetic  operation  and  assigning  it  a 
normalized  time  based  on  average  Assembler  instruction  times  for  the 
system.  This  method  assumes  the  following: 

(1)  the  relative  times  between  different  Assembler  and  Fortran 
instructions  are  about  the  same  (i.e.  a division  takes  about 
5.8  times  as  long  as  an  addition  regardless  of  the  code 
used) ; and 

(2)  the  time  involved  in  transfer  of  control,  variable 
initializations  and  other  overhead  are  small  compared 
to  the  time  required  to  perform  arithmetic  operations. 

f 
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STEP 


CAM)  LISTING 


1 //**  (Resource  Card  if  Needed) 

(Charge  Card) 

2 //_EXEC_FORTGC 

3 //FORT.  SYSLINJDD_SPACE  - (CYL,  (2,1)) 

4 (FIRST  CARD  IN  FORTRAN  DECK) 


. CALL  STIMER  (initialize  clock) 

CALL  SUBROUTINE  ABSTD  (to  be  timed) 
CALL  TTIMER  (ITIME)  (stops  clock) 
TIME  - ITIME*26. 04116  (in  usee) 
WRITE  (6, XX)  TIME 
XX  FORMAT  (output) 

(LAST  CARD  IN  FORTRAN  DECK) 

5 //* 

6 //__EXEC_ASMFCLG,  PAKM.  ASM  - 'NO  DECK, 

LOAD',  PARM.LKED  - 'LIST,  LET, 
XREF,  CALL' 

7 STIMER  CSECT 

BASRG  12 

STIMER  TASK,  TUINTVL  - BIG 
BYE 

ENTRY  TTIMER 
TTIMER  BASRG  12 

L 2,  0(1) 

TTIMER 

L 1,  BIG 

SR  1,  0 

ST  1,  0(2) 

BYE 

BIG  DC  0F'0',  X'0FFFFFFF' 

END 

8 / /LKED . SYSLMOD  DD  SPACE  - (CYL,  (5,5,1)) 


9 / /LKED . SYSUT I_DD  SPACE  - (CYL,  (5,5,1)) 


STEP 


CARD  LISTING 


10  / /LKED . SYSLIB_DD_D5N  - SYSl.FORTLIB,  DISP-SHR 
//  DD  DSN  - UTEP,  FORTLIB,  DISP-SHR 

//  DD  DSN  - UTEP,  CALCOHP . LOAD , DISP-SHR 

11  //GO.FT06F001  DD  SYSOUT  - A 

12  //GO.FT05F001  DD  * 

13  (FORTRAN  DATA  CARDS  IF  NEEDED) 

14  II* 

15  (FORTRAN  JCL  CARDS  IF  NEEDED) 

16  // 


Figure  5.5-1  (Cont.)  TIMING  ROUTINE 


Table  5.5-1  la  a list  of  typical  arithmetic  operations  and  the 
average  Assembler  Instruction  time  In  psec  for  the  IBM  360/65  system 
with  H-level  compiler.  Instruction  times  are  based  on  double-precision 
real  variable  operations.  Column  4 lists  the  times  normalized  to 
that  required  by  addition  (TA) . Times  which  had  to  be  estimated  are 
indicated . 

In  summary,  three  types  of  timing  schemes  can  be  used  to  indicate 
which  routine  is  computationally  faster  than  another.  The  advantages 
and  disadvantages  of  each  method  are  presented  in  Fig.  5.5-2.  All 
three  methods  were  used  in  timing  the  routines  for  computing  KgTD  by 
Eq.  5.1-1  and  Kj^  by  Eq.  5.1-2.  Results  of  the  time  analyses  are 
presented  in  Chapter  6. 


OPERATION 

SYMBOL 

TIME (psec) 

NORMALIZED  TIME 

Addition 

+ 

2.45 

1ta 

Subtraction 

- 

2.45 

1ta 

Division 

T 

14.10 

5.8  Ta 

Multiplication 

* 

7.60 

3-1ta 

Exponentiation 

DEXP 

145.00 

59.2  Ta 

Square  Root 

DSQRT 

5.85 

2.4  Ta 

Integer  Power 

**INT 

(7. 6)  INT 

(3.1) (INT)  T. 

A 

Real  Power 

**REAL 

145.00 

59.2  Ta 

Logical  Branch 

IF 

5.85 

2.4  Ta 

Loop 

DO 

92 

37.6  Ta 

Note:  (1)  Symbol  INT  represents  the  integer  power  to  which 

a real  # is  raised. 

(2)  Ta  is  the  time  for  1 addition  in  psec.  - i.45  psec 
for  Assembler. 
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1 


GO- STEP  TIMER 


Advantage  - Good  overall  accuracy,  easy  to  use, 
directly  related  to  computer  usage 
costs. 

Disadvantage  - Only  measures  complete  program, 
fluctuates . 


CLOCK-TYPE  TIMER 


Advantage  - Good  accuracy,  can  time  selected 

portion(s)  of  routine,  can  be  used 
more  than  once  in  a given  run. 

Disadvantage  - Requires  some  programming  knowledge, 
fluctuates  from  run  to  run. 


ASSEMBLY-COMPARISON 


Advantage  - Good  relative  measure,  requires  no 

programming,  measures  only  arithmetic 
operations . 


Disadvantage  - Not  accurate  in  computer  time  sense, 
requires  estimates , excludes  overhead 
times . 


I 


Figure  5.5-2  TIME  ANALYSIS  METHODS 
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5.6-lc 


C3)  Percent  Absolute  Relative  Error  - 
E%  « 10Q  Er- 

The  quantity  represents  the  standard  and  represents  the  appraxi- 
mation.  Absorption  coefficient  K,  and  transmittance  t,  calculations 
computed  in  the  standard  fashion  are  subscripted  "STD";  those  done  by 
the  approximate  method  are  subscripted  "MOD." 

Absolute  error,  E^,  is  sometimes  called  a "deviation."  It 
measures  the  difference  between  two  quantities.  Because  transmittance 
is  always  between  0 and  1,  the  order  of  magnitude  of  E^  will  also  in- 
dicate the  number  of  decimal  digits  of  accuracy  obtained  with 
Relative  error,  Er>  measures  the  fractional  part  an  estimate  calcula- 
tion is  "off"  from  the  standard.  Percent  absolute  relative  error, 

E„,  is  simply  stating  E is  terms  of  a percentage. 
a r 

To  illustrate  the  concept  of  E^  and  Ey  as  they  could  apply  to 
transmittance  calculations  consider  the  two  cases  below: 

(1)  Case  I Let 

TSTD  and  TMOD  " *99985. 

Therefore, 

EA  - 5 x 10“5, 

and 

E - 5.0005  x 10-5. 
r 
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>*  .Q001  and  t. 


. 000Q5 


TSTD 


MOD 


Again 


E - 5 x 10 

A 


but 

E - 5 x 10_1, 
r 

which  is  an  increase  of  4 orders  of  magnitude  in  absolute  relative 

error,  even  though  absolute  error  has  remained  the  same. 

Conversely,  if  the  desire  had  been  to  keep  relative  error  the 

same  as  in  Case  I,  the  absolute  error  in  Case  II  would  have  to  be 
_9 

on  the  order  of  10  which  is  an  increase  of  4 orders  of  magnitude. 

Caution  must  be  exercised  in  error  predictions  for  the  reasons 
illustrated  above.  Error  must  be  evaluated  within  the  total  context 
in  which  calculations  are  to  be  used.  Because  the  ultimate  use  of 
absorption  coefficient  calculations  are  to  compute  transmittances,  it 
will  be  within  this  scope  that  error  for  the  approximation  will  be 
analysed . 

Beer’s  Law  (Chapter  2)  states  that,  for  monochromatic  radiation, 
transmittance  is 


t - EXP  [-KU], 


5.6-2 


where  U is  the  amount  of  gas  that  the  radiation  must  pass  through 
Hence,  the  absolute  error  is 


5.6-3 


ea- 


tstd  tmod 


EXPI-KStdUJ  - EXP I-V] 


KSTD  an^  ^HOD  are  re*ate<*  by  their  absolute  error: 
CA  “ |KSTD  ~ ^MOD  * 


5.6-4 


where  the  symbol  will  distinguish  absorption  coefficient  error  from 
transmittance  error,  E^.  Combining  Eqs.  5.6-3  and. 5.6-4,  can  be 
expressed  in  terms  of  Tg^,  U,  and  as  follows: 


STD 


1 - EXP£eAU] 


5.6-5 


Note  that 


1 - EXPl+eAU] 


E . 
r 


5.6-6 


It  can  be  seen  from  Eq.  5.6-5  that  EA  will  vary  depending  on  the  value 
of  as  well  as  the  quantity  (eAU) . 

If  three  digits  of  accuracy  (rounded)  are  desired  in  transmittance 
then  a bound  can  be  placed  on  the  value  of  EA: 

E.  < 5 x 10  4. 

A — 

Then  Eq.  5.6-6  can  be  written  to  reflect  this  upper  limit  as 


5 x 10~4 
> 

tstd 


1 - EXPl+eAU] 


5.6-7 
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It  can  be  seen  that  the  inequality  will  be  most  stringent  for  trans- 
mittance values  close  to  unity  and  moat  relaxed  for  values  near  zero 


Table  5.6-1  is  a listing  of  the  limits  placed  on  the  quantity  (e^U) 
for  various  values  of  TgTD.  Computed  values  of  (eAU)  for  the  4.3ym 
band  are  given  in  Chapter  6. 


5 x 10 
TSTD 


-4 


1 - EXP[+e^U] 


Table  5.6-1  Bounds  on  |e.U|  for 


E.  - 5 x 10 
A 


A 

-4 
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CHAPTER  VI 


RESULTS  AND  APPLICATIONS 


6.1  Preface 

The  various  methods  of  quantifying  computation  speed  and 
accuracy  for  the  approximation  to  the  Lorentz  line  absorption 
coefficient  were  outlined  in  Chapter  5.  Sections  6.2  and  6.3  of 
this  chapter  present  the  time  and  accuracy  measurements  obtained 
by  testing  the  subroutines  "ABSTD"  and  "ABSMOD"  listed  in  Appendix 
B.  All  test  runs  were  conducted  on  the  University  of  Texas  at 
El  Paso's  IBM  360/65  computer.  Results  using  Chebyshev  economiza- 
tion to  replace  the  Taylor  series  expansion  (discussed  in  Section 
3.4)  is  also  contained  in  this  chapter.  Finally,  a brief  section 
is  included  on  the  expected  results  of  applying  the  approximation 
to  a transmittance  "band"  problem. 

6.2  Results  of  the  Time  Analysis 

To  calculate  the  time  involved  in  executing  subroutines  "ABSTD" 
and  "ABSMOD"  (Appendix  B)  for  monochromatic  radiation,  each  timing 
method  discussed  in  Section  5.5  was  used.  In  all  cases  a signifi- 
cant time  savings  was  measured  when  the  absorption  coefficient  due 
to  Lorentz  broadening  was  calculated  by  the  approximation  express  in 
Eq.  5.1-2. 

All  timing  methods  sought  to  measure  only  the  calculation  time 
involved  in  computing  K.  Every  effort  was  made  to  subtract  out  the 
time  involved  in  data  retrieval,  routine  initializations,  and  other 
programming  overheads. 


The  "GO-STEP"  computation  times  were  obtained  by  making 
computer  runs  with  and  without  different  sections  of  the  absorption 
coefficient  program  and  subtracting  the  CPU  "GO-STEP"  times.  The 
types  of  runs  timed  were: 

(1)  Program  without  subroutines  "ABSTD"  and  "ABSMOD"; 

(2)  Program  without  subroutines  "ABSTD"  and  "ABSMOD"  and 
also  without  the  "DATA  REORGANIZATION"  portion  of 
the  main  routine; 

(3)  Program  containing  only  subroutine  "ABSTD"; 

(4)  Program  using  only  subroutine  "ABSTD"  modified  to 
return  control  as  soon  as  line  parameters  were  read 
into  core;  and  finally, 


(5)  Program  containing  only  subroutine  "ABSMOD" 


If  T^  - T,.  are  defined  as  the  "GO-STEP"  times  for  (1)  - (5)  above 

and  T_„_  and  T„„.  represent  the  time  required  to  calculate  K--..  and 
bTD  MUD  bJLD 


respectively,  then 


TSTD  T3  “ T2" 


TMOD  T5  “ T2' 


The  time  involved  in  reorganizing  the  line  parameter  data  for  use 
with  subroutine  "ABSMOD"  is  given  by 


T = T — T 

DATA  1 2’ 


Using  the  above  definitions,  timing  runs  were  made  using  a 


standard  33  — level  atmosphere  (see  Table  6.3-1)  for  channels  11—15 
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25 

[Si8sila  ].  Representative  times  are  given  in  Table  6.2-1.  The 
times  listed  for  TgTD  and  TMQD  represent  the  best  estimates  of  CPU 
time  required  to  compute  KgTD  and  respectively. 


CHANNEL 

NUMBER 

LINE 

DENSITY 

T 

STD 

(sec) 

T 

MOD 

(sec) 

T 

DATA 

(sec) 

T 

STD 

T 

MOD 

11 

993 

15.38 

1.59 

5.18 

9.67 

12 

1058 

16.04 

1.72 

6.02 

9.32 

13 

921 

12.82 

1.44 

5.10 

8.90 

14 

1597 

44.14 

2.83 

10.68 

15.59 

15 

504 

5.15 

1.01 

3.42 

5.09 

NOTE:  Subroutine  "ABSMOD"  used 
and  Taylor  Series . 

7 terms 

in  both  Binomial 

Table  6.2-1  "GO-STEP"  TIME  C/ART 

Note  that  the  numbers  of  spectral  lines  which  contribute  to  the 
absorption  are  listed  also. 

For  channels  11,  12,  and  13  which  all  contain  around  1000 
spectral  lines,  was  calculated  about  9 times  faster  than 

As  the  number  of  lines  in  the  channel  increase  (decrease)  the  ratio 
increases  (decreases).  A ratio  of  over  15  was  achieved  in  channel 
1A  where  the  number  of  lines  exceed  1500. 

If  the  user  is  willing  to  make  a small  study  of  the  line 
parameters  contained  in  the  spectral  band  of  interest  and  compute 
the  "DATA  REORGANIZATION"  portion  of  the  routine  separately,  a 
significant  savings  in  computer  time  can  be  realized  in  calculating 
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K.  A lesser  sayings  will  be  attained  if  the  data  is  reorganized  in 
the  same  routine  that  K^0D  is  calculated.  For  monochromatic  radiation. 
It  is  always  possible  to  reorganize  the  line  parameter  data  Separately. 
Unfortunately , for  polychromatic  considerations , only  a portion  of 
the  data  can  be  reorganized  this  way.  Section  6.4  will  address  this 
case. 

Measurements  using  the  "GO-STEP"  timer  are  difficult  to  adapt 
to  a particular  portion  of  a computer  program  (i.e.  subroutines 
"ABSTD"  and  "ABSMOD")*  This  method  requires  numerous  program  executions 
with  and  without  selected  sections  which  are  to  be  timed.  Because  of 
this,  the  "GO-STEP"  method  is  awkward  and  inconvenient  to  implement, 
and  hence,  not  convincingly  accurate. 

In  an  effort  to  measure  the  execution  times  of  subroutines  "ABSMOD" 
and  "ABSTD"  more  precisely,  the  "CLOCK"  timer  (see  Fig.  5.5-1)  was  also 
used.  The  execution  time  of  each  subroutine  was  measured  independent- 
ly from  the  main  program  because  the  "CLOCK"  timer  enables  the  user 
to  start  and  stop  the  clock  whenever  desired.  Using  the  "CLOCK"  timer 
routine  presented  in  Chapter  5,  similar  timing  runs  were  made  to  those 
in  Table  6.2-1.  The  "CLOCK"  times  for  each  subroutine  and  the  "DATA 
REORGANIZATION"  section  are  given  in  Table  6.2-2.  TM0D  was  found  to 
be  about  14  times  faster  than  TgTD  for  1000  spectral  lines  considered. 

As  the  number  of  lines  approached  1600  the  ratio  of  computation  times 
(column  6,  Table  6.6-2)  exceeded  16. 
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CHANNEL 

NUMBER 

LINE 

DENSITY 

T 

STD 

(sec.) 

T 

MOD 

(sec.) 

T 

DATA 

(sec.) 

T 

STD 

T 

MOD 

11 

993 

16.15 

1.08 

5.53 

14.95 

12 

1058 

16.35 

1.13 

6.13 

14.46 

13 

921 

16.04 

1.10 

5.06 

14.58 

14 

1597 

30.12 

1.88 

9.86 

16.02 

15 

504 

5.51 

0.62 

3.35 

8.89 

NOTE : 

Subroutine  "ABSMOD"  used  7 terms  in  both 
Binomial  and  Taylor  Series. 

TABLE  6.2-2  "CLOCK"  TIME  CHART 

Tables  6.2-1  and  6.2-2  suggest  that  TgTD  and  TM0D  are  strongly 
dependent  on  the  number  of  spectral  lines  present  within  the  band 
limits.  This  dependence  can  be  quantized  using  the  method  of 
"ASSEMBLER"  time  comparison  which  was  outlined  in  section  5.5.  If 
it  is  assumed  that  0.5%  of  the  spectral  lines  will  be: 

(\>  - v ) < 0.0525 
o — 

wave  numbers , and  only  7 terms  of  the  Binomial  and  Taylor  Series  are 
used,  then  for  the  program  in  Appendix  B, 

TSTD  = t 38204 -1  + 10490.86  I]  Ta 

and 


- 11325982.2  + 23.8978  I]  T , 

MOD  A 


Where  I ia  the  total  numher  of  spectral  lines  considered,  and  is 
the  time  in  p sec  to  complete  one  assembler  addition  instruction.  The 
ratio  of  to  ^jjqD  f°r  ^ values  of  500,  1000,  and  1500  are  listed 
in  Table  6.2-3.  Values  for  TgTD  and  T^  ^ are  given  as  functions  of 
Ta.  Tjjqd  *8  predicted  to  be  7.8  times  faster  than  T^p  for  the  1000 
line  case. 


I 


I 

VALUE 

* T 

STD 

(sec.) 

T 

MOD 

(sec.) 

tstd 

tmod 

500 

5.28  T. 

A 

1.34  T. 

A 

3.94 

1000 

10.53  T. 

A 

1.35  Ta 

A 

7.80 

1500 

15.77  T. 

A 

1.36  Ta 

A 

11.59 

TABLE  6.2-3  "ASSEMBLER"  TIME  CHART 


6.3  Results  of  Accuracy  Analysis 


In  section  5.6  various  definitions  of  error  were  defined.  Table 


5.6-1  lists  typical  bounds  on  the  absolute  error  in  K,  e^,  times  the 
gas  amount,  U.  These  bounds  assure  3 place  accuracy  in  transmittance 
(vhen  rounded  to  the  nearest  ten-thousandths) . The  inequality  which 


must  be  satisfied  is  rewritten  here  for  convenience  as 


5.0  x 10 
TSTD 


1 - EXP  H eA  u] 


6. 3—1 
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where  x 


la  the  transmittance  through  gas  amount  U calculated  using 


STD 


the  standard  Lorentz  line  absorption  coefficient,  and  is  the  abso- 
lute error  in  absorption  coefficients  as  defined  by  Eq.  5.6-4. 

To  calculate  K,  U,  and  e^,  a standard  atmospheric  profile  contain- 
ing 33  levels  was  used.  The  profile  which  represents  typical  tempera- 
tures and  pressures  for  a midiattitude  - summer  climate  is  listed  in 
Table  6.3-1.  Altitude  ranges  from  sea-level  to  100  KM.  Pressures 
are  listed  in  both  millibars  and  atmospheres. 

Gas  amounts  were  calculated  for  every  layer  between  each  level 
of  the  profile.  All  seven  gas  amounts  were  computed,  however,  Table 
6.3-2  only  lists  the  most  important  layers  for  the  four  gases  found 
in  the  4.3  pm  band  0^0 , N^O,  CO).  It  should  be  noted  that  for 

all  gases  listed,  the  maximum  gas  amount  is  found  in  layer  //I  (0  to  1 
km) . The  only  exception  in  all  seven  gases  analyzed  was  ozone  which 
had  its  maximum  concentration  occur  between  21  to  22  km  (layer  22) . 
Column  2 of  Table  6.3-2  lists  the  lower  and  upper  altitude  limits  for 
each  layer  given. 

Using  the  program  in  Appendix  B values  for  KgTD  and  were 

computed  for  all  channels  in  the  4.3  pm  band.  Using  7 terms  in  the 
Binomial  and  Taylor  Series  expansions  for  computation  of  K^QD,  three 
place  accuracy  in  transmittance  was  more  than  assured  in  all  cases. 

To  insure  this  result  the  following  analysis  was  performed. 

Let 


in 


5 x 10 


-4 


STD 


6.3-2 
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(KM) 

PRESSURE 

(mb)  (ATM) 

TEMPERATURE 

(*K) 

Q 

1013.0 

0.999 

294 

1 

902.0 

0.890 

290 

2 

602.0 

0.791 

285 

3 

710.0 

0.700 

279 

4 

628.0 

0.619 

273 

5 

554.0 

0.546 

267 

6 

487.0 

0.480 

261 

7 

426.0 

0.420 

255 

8 

372.0 

0.367 

248 

9 

324.0 

0.319 

242 

10 

281.0 

0.277 

235 

11 

243.0 

0.239 

229 

12 

209.0 

0.206 

222 

13 

179.0 

0.176 

216 

14 

153.0 

0.151 

216 

15 

130.0 

0.128 

216 

16 

111.0 

0.109 

216 

17 

95.0 

0.093 

216 

18 

81.2 

0.080 

216 

19 

69.5 

0.068 

217 

20 

59.5 

0.058 

218 

21 

51.0 

0.050 

219 

22 

43.7 

0.043 

220 

23 

37.6 

0.037 

222 

24 

32.2 

0.031 

223 

25 

27.7 

0.027 

224 

30 

13.2 

0.013 

234 

35 

6.52 

0.643  x 10~2 

245 

40 

3.33 

0.328  x 10"2 

258 

45 

1.76 

0.173  x 10'2 

270 

50 

0.951 

0.938  x 10'3 

276 

70 

0.061 

0.662  x 10“4 

218 

100 

0.0003 

0.296  x 10-6 

210 

TABLE  6.3-1  STANDARD  33-LEVEL  ATMOSPHERE  [Adapted  from  La  Rocca16] 
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Layer 

No. 

Alt 

AZ  (KM) 

H2° 

C02 

n2o 

CO 

1 

| 

0-1 

1.152 

28.697 

.024 

.006 

2 

1-2 

.747 

25.853 

.022 

.005 

3 

2-3 

.458 

23.785 

.020 

.005 

4 

3-4 

.255 

21.200 

.017 

.004 

5 

4-5 

.141 

19.131 

.016 

.004 

6 

5-6 

.079 

17.322 

.014 

.004 

7 

6-7 

.048 

15.770 

.013 

.003 

8 

7-8 

.028 

13.961 

.011 

.003 

9 

8-9 

.016 

12.409 

.010 

.002 

10 

9-10 

.009 

11.117 

.009 

.002 

15 

14-15 

8.9  x 10"5 

5.946 

.005 

.001 

20 

19-20 

4.6  x 10 

2.585 

.002 

5.9  x 10"4 

25 

24-25 

6.3  x 10"5 

1.163 

9.8  x 10~4 

2.6  x 10-4 

30 

45-50 

7.1  x 10“6 

.209 

1.7  x 10~4 

4.7  x 10"5 

31 

50-70 

3.9  x 10-6 

.228 

1.9  x 10~4 

5.2  x 10~5 

32 

70-100 

1.3  x 10'7 

.017 

1.4  x 10~5 

3.9  x 10-6 

Note:  Units  of  gas  amounts  are  (ATM-CM) . 

TABLE  6.3-2  COMPUTED  GAS  AMOUNTS 
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which  is  the  most  restrictive  case  of  Eq.  6.3-1.  In  all  cases  the 

maximum  yalue  of  e.U  was  obtained  between  0 and  1 km.  Values  for 

A 


Jin 


5 x 10 
TSTD 


are  given  in  Table  6.3-3  to  6.3-6.  The  transmittance,  TgTD,  was  com- 
puted for  the  same  layer  in  which  the  gas  amount  was  found  to  be  max- 
imum. 

Even  though  the  error  tolerances  in  transmittance  were  not  breached, 
Table  6.3-4  shows  a definite  decrease  in  accuracy  for  channels  13-15 
(CX^).  This  is  caused  by  the  Taylor  series  expansion  of  EXP[X]  and  the 
variations  of  the  magnitude  of  X.  The  use  of  Chebyshev  economization 
will  correct  the  skewed  error  distribution. 

Recall  the  Taylor  series  used  in  the  approximation  (Eq.  3.4-6): 


EXPlX] 


EXP 


- [ 
5T*  L 


pvv(vv  j 
> <vv  1 ° 

kT  J 


<Ea-El 


n! 


n"0 


6.6-3 


where  the  bracketed  term,  X,  is  approximated  using  (N+l)  terms.  When 
X becomes  large,  more  terms  are  needed  to  approximate  EXP[X]  to  keep 
the  same  degree  of  accuracy.  To  keep  the  magnitude  of  X small,  strict 
restrictions  have  to  be  placed  on  constants  E and  W . 

a a 

The  value  of  X is  determined  by  two  primary  factors: 
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Churns  1 
No. 

[6aV1mL 

ERROR 

BOUND 

tstd 

11 

3.439  x 10"10 

4.999  x 10"4 

.9999929 

12 

3.228  x 10'10 

4.999  x 10-4 

.9999690 

13 

1.013  x 10'9 

4.999  x 10-4 

.9999750 

14 

7.756  x 10-5 

4.999  x 10-4 

.9999224 

15 

3.067  x 10"10 

4.999  x 10-4 

.9999926 

*U  ■ 1, 

MAX 

,152  atm-cm  of  I^O 

I 


TABLE  6.3-3  ERROR  ANALYSIS  - WATER  VAPOR 
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TABLE  6.3-4  ERROR  ANALYSIS  - CARBON  DIOXIDE 


CHANNEL 

NO. 

L * JMAX 

ERROR 

BOUND 

tstd 

11 

2.323  x 10-8 

5.000  x 10'4 

.9996881 

12 

2.547  x lO-7 

5.095  x 10-4 

.9810430 

13 

2.130  x 10“4 

1.225  x lO-3 

.4078124 

14 

1.935  x 10"3 

2.509  x 10-1 

.0019920 

15 

3.600  x 10-1 

->  00 

0.0  x 10'" 

*11.. _ “ 28.697  atm  -cm  of  C0o 

MAX  x 

T»bl«  6.3-4  ERROR  ANALYSIS-CARBON  DIOXIDE 


CHANNEL 

NO. 

[eA  vr 

U L JHAX 

ERROR 

BOUND 

tstd 

11 

8.522  x 10~6 

5.354  x 10~4 

.9336696 

12 

6.295  x 10“5 

6.746  x 10“4 

.7409964 

13 

6.514  x 10~5 

7.255  x 10“4 

. 6890441 

14 

1.413  x 10~9 

4.999  x 10-4 

.9999888 

15 

1.865  x 10~9 

4.999  x 10“4 

.9999779 

024  atm- cm  of  N^O 


CHANNEL 

NO. 

MAX 

ERROR 

BOUND 

tstd 

11 

0.0 

6.58  x 10-4 

.8117192 

12 

7.263  x 10-9 

5.001  x 10~4 

.9995790 

13 

6.617  x 10“10 

.4999  x 10-4 

.9999809 

14 

1.242  x 10"15 

4.999  x 10~4 

1.0 

15 

N/A 

N/A 

N/A 

*UW , ^ - .006  atm-cm  of  CO 
MAX 


Table  6.3-6  ERROR  ANALYSIS-CARBON  MONOXIDE 
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X - AXjXj, 


where  A is  a constant,  and 


Xl-  lEa-Ei'* 


x2  - |wx  - vg. 


Values  of  range  from  a low  of  0 to  a maximum  in  excess  or  4700 
wavenumbers.  If  E&  is  picked  to  minimize  X^,  variations  may  still 
exceed  2000  wavenumbers.  Based  on  a standard  temperature  profile, 

X2  will  vary  from  -0.3  to  +0.3.  This  implies  that  X,  using  appro- 
priate conversion  factors  in  A,  will  fall  between  +4.5  depending  on 
energy  and  temperature  values.  Values  of  X which  reach  this  magnitude 
may  require  up  to  20  terms  in  the  Taylor  expansion  and  this  undermines 
the  inherent  quickness  of  the  approximation. 

A study  of  the  line  parameter  data  revealed  the  fact  that  high 
energy  values  are  associated  with  low  line  intensity  values.  Hence 
it  is  possible  to  limit  the  value  of  X,  by  placing  a maximum  on  the 
ground-state  energy.  Analysis  indicates  that  7 decimal  place  accu- 
racy in  K„__.  is  maintained  when  the  maximum  energy  value  is  limited 
to  3250  wavenumbers.  Using  this  maximum,  X can  be  reduced  to  values 
in  the  range  of  -2.5  to  +2.5.  Using  Chebyshev  coefficients  to 
equally  distribute  the  error,  a greater  stability  in  the  approxima- 
tion is  attained. 

Table  6.3-7  lists  the  absolute  error  in  absorption  coefficients, 

j^HOD  “ KSTd|*  where  Kjjqd  *8  calculated  by  both  Chebyshev  and  Taylor 
expansions.  Note  that  the  Chebyshev  error  in  K,  expressed  in  atmos- 
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pharlc  centimeters , is  as  good  or  better  than  the  Taylor  error.  The 


r 


* 


\ 


bast  correction  occurs  In  channel  15  for  CO^. 

In  Table  6.3-8,  Chebyshev  coefficients  for  the  EXPlX]  expansion 
are  given  for  reference.  Coefficients  are  listed  for  various  ranges 
in  X and  for  different  numbers  of  terms.  The  standard  form  of  the 
Chebyshev  expansion  Is 

expix]  - cQ  + c_jX  + c/  + c3  x3  + ....  + cnxn, 

where  (n  + 1)  terms  are  used. 


6 . 4 Application  to  the  Band  Problem 

The  inherent  importance  of  a fast  and  accurate  approximation  to 
the  Lorentz  line  absorption  coefficient  is  illustrated  more  in  the 
light  of  transmittance  calculations.  Thus  far  transmittance  has  been 
discussed  in  terms  of  only  monochromatic  radiation  which,  in  reality, 
can  not  actually  be  measured.  Because  present  instrumentation  is 
restricted  to  opeating  in  a spectral  band  regardless  of  the  resolution 
achieved,  transmittance  is  often  desired  over  a range  of  frequencies. 

The  polychromatic  transmittance  is  normally  expressed  as  an  aver- 
age of  the  monochromatic  values  taken  over  a band.  Mathematically, 


T 


1_ 

V 

o 


t (v,T ,P)dv, 

V 

o 


6.4-1 


where  v is  the  spectral  band  width  and  t(v,T,P)  is  the  monochromatic 
o 

transmittance.  Regardless  of  the  transmittance  desired,  the  Lorentz 
line  absorption  coefficient  can  be  calculated  by  the  approximation 
defined  in  Eq.  3.3-8.  Accuracy  of  x will  be  similar  to  that  of  t. 
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Channel /Gas 
11/HjO 

12/HjO 

13/H20 

14/HjC 

15/H2° 

n/co2 

12/COj 

13/C02 

14/C02 

15/C02 

ii/n2o 

12/N20 

13/N20 

14/N20 

15/N20 

11/CO 
12  /CO 
13/CO 
14/CO 


Chebyehev  Error*  Taylor  Error 


2.2  x 

3.0  x 
8.9  x 

1.3  x 

3.5  x 

1.1  x 

9.0  x 

7.0  x 

6.6  x 
8.9  x 

1.3  x 

2.6  x 

3.2  x 

5.8  x 

5.7  x 

1.9  x 

1.2  x 

7.1  x 
2.0  x 


10 


-11 


10 


-10 


10 


-12 


10 


-9 


10 


-12 


10 


-10 


10 


-9 


10 


-8 


10 


-5 


10 


-5 


10 


-5 


10 


-3 


10 


-4 


10 


-8 


10 


-9 


10 


-7 


10 


-6 


10 


-9 


10 


-13 


2.9  x 10 


-10 


3.2  x 10 


-10 


8.6  x 10 


-10 


6.7  x 10 


6.3  x 10 


r5 

-6 


8.0  x 10 


-10 


8.8  x 10 


-9 


7.4  x 10 


-6 


6.7  x 10 


-5 


1.2  x 10' 


-2 


3.4  x 10 

2.6  x 10 

2.7  x 10 

5.8  x 10 


-4 


-3 


-3 


-8 


-8 


7.7  x 10 


2.6  x 10-6 


1.2  x 10 


-6 


1.1  x 10 


-7 


2.0  x 10 


-13 


Note:  Chebyshev  expansion  used  9 terms  for  X on  the  Interval 
1-2.7  to  +2.7] 


Table  6.3-8  CHEBYSHEV  COEFFICIENTS  FOR  EXPlX] 
[Adapted  from  Gibson^] 
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Tha  method  of  calculation,  and  hanca,  tha  apaad  of  tha  approximation, 
will  chanfa  howavar. 

Za  tha  monochromatic  caaa,  tha  approximation  achlavaa  ita  graat 
apaad  advantaga  by  alloying  thoaa  quantities  dapandant  on  llna  para- 
aatara  vQi,  a^,  SQ^  and  to  ba  calculatad  aaparata  from  K.  For 
a band  transmittance,  on  tha  otharhand,  quantit iaa  dapandant  on 
can  not  ba  calculatad  prior  to  tha  calculation  of  K bacauaa  of  tha 
fraquancy  atapping  that  muat  ba  dona. 

For  axampla,  tha  monochromatic  "DATA  REORGANIZATION"  routlna 
computaa 


W- 


vhara  I la  tha  nuabar  of  linaa.  Tha  raaulta  ara  atorad  in  a (NxM) 
array  for  aach  gaa  and  naad  ba  computad  only  onca  for  a channal.  For 
tha  polychromatic  Inatanca,  tha  tarn 

fm(Tt)  a (Y±2  ♦ aV(B+1)  , 


vhara 


v - v 


Yi* 


oi 


oi 


2 

and  a ia  daflnad  in  Chaptar  3,  can  not  ba  computad  onoa  and  forgottaa. 
Bacauaa  tha  valuaa  of  v ara  atappad  acroaa  tha  apactral  band  In  small 
incramanta,  fffl(Yj)  vill  changa  in  valua  aach  tlma  a nav  \>  ia  uaad. 
Oibaon  1 9]  haa  aatlmatad  tha  pradictad  TSTJ)  and  timaa  for 
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the  polychromatic  case.  The  estimates  are  baaed  on  the  "ASSEMBLER" 
method  explained  previously  and  projected  calculations  which  will  be 
required.  At  the  time  of  this  writing  the  approximation  has  not 
actually  been  applied  to  the  polychromatic  problem  and  verified  by 
computational  analysis. 

Letting  P represent  the  number  of  frequency  increments  taken 
within  a band  of  width  vq,  the  estimated  times  are  given  by  Gibson  as 

T^  * (113  PI  + 1245L  + 415  PL)TA 

and 

Tstd  A (15  PIL  + 906  L + 62  IL)TA> 

where  there  are  "I"  lines  considered  within  the  band  and  "L"  is  the 
number  of  stratified  levels  of  the  atmosphere.  TA  is  the  time  required 
to  perform  an  assembler  addition  instruction. 

Using  these  predictors  for  TMQD  and  TgTD»  Table  6.4-1  was  developed. 
The  frequency  increment  (P)  used  was  based  on  a band  width  of  25  wave- 
numbers  and  a step  size  of  0.02  cm  \ 


F 
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I 

value 

T 

STD 

(sec) 

tmod 

(sec) 

tstd 

,£mod 

500 

310.43 

87.78 

3.53 

1000 

620.83 

158.40 

3.92 

1500 

931.22 

229.03 

4.06 

2000 

1241.62 

299.65 

4.14 

P - 1250 
L - 33 


T _ - (113  PI  + 1245  L + 415  PL)T. 

MUD  A 

tstd  * (15  PIL  + 906  L + 62  il>ta 


Table  6.4-1  POLYCHROMATIC  TIME  CHART 


CHAPTER  yil 


CONCLUDING  REMARKS 

7.1  Summary 

Attenuation  of  radiant  energy  by  the  earth's  atmosphere  can  be 
described  mathematically  by  the  Lambert-Beer  Lav  (see  Sec.  2.4).  In 
the  Infrared  spectrum,  the  attenuation  Is  due  to  molecular  absorption 
by  water  vapor,  carbon  dioxide,  ozone,  nitrous  oxide,  methane,  carbon 
monoxide,  and  sometimes,  oxygen.  The  absorption  coefficient,  K,  Is 
composed  of  a strength  and  a line  shape.  The  line  shape  In  the  lover 
regions  of  the  atmosphere  Is  primarily  due  to  pressure  broadening. 

ThA  standard  formula  for  the  absorption  coefficient  Is  glvsn  by  Eq. 
2.3-7a. 

The  Lorentz  line  absorption  coefficient  can  be  written  In  an 
alternate  form  as  given  by  Eq.  3.3-8  and  approximated  using  Eq.  5.1-2. 
Using  this  approximation,  K can  be  computed  about  9 times  faster  than 
the  standard  method,  while  retaining  at  least  three  digits  of  accuracy 
in  corresponding  transmittance  calculations.  Accuracy  can  be  improved 
even  more  by  replacing  the  Taylor  series  coefficients  with  those  of 
Chebyshev.  The  approximation's  development  and  its  testing  are  pre- 
sented In  Chapters  3,  5,  and  6. 

Chapters  2 and  4 are  intended  to  lend  emphasis  to  the  problem 
of  atmospheric  transmittance  as  a whole.  Often  the  theoretical 
developer  neglects  the  mundane  problems  of  atmospheric  units  and 
data.  Hopefully,  these  chapters  will  enable  the  novice  atmospheric 
scientist  to  better  understand  the  complexities  with  which  he  Is  dealing. 


7.2  Use  of  the  Approximation 

The  approximation  to  the  Lorentz  line  absorption  coefficient 
will  sufficiently  retain  the  high  resolution  calculations  of  the 
standard  line-by-line  method  and  simultaneously,  greatly  reduce  the 
computational  time  required.  Continued  investigation  is  needed  into 
applying  the  approximation  to  a polychromatic  transmittance  problem. 
Prom  results  obtained  thus  far,  the  adversion  of  using  line-by-line 
methods  to  calculate  "band"  transmittances  may  be  eliminated. 
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APPENDIX  A 


DATA  ROUTINE  FOR  USE  WITH  AFCRL’S 
LINE  PARAMETER  COMPILATION 


FORTRAN  IV  0 LEVEL  21 


MAIN  DATE  » 77108  ll/31/i 


0001 

0002 

0003 

0004 

0005 

0006 
0007 
00  08 

0009 

0010 
0011 
0012 

0013 

0014 

0015 


0016 

0017 


c 

c 

c 

c 

r. 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


c 

c 

c 


4.3  MICRON  BAND 
MICHAEL  T.  POTTER 


THIS  PROGRAM  USES  SUBROUTINE  "CHPEAD"  AND  "MCCL ATCHEV • S"  TAPE 
TO  OUTPUT  SPECTRAL  LINE  DATA  WITHIN  SELECTED  8ANOWIOTHS. 

IN  ADDITION,  THIS  PROGRAM  WILL  ALSO  B'JILO  PERMANENT  DATA  FILES. 
PRESENT  PROGRAM  INITIALIZATIONS  ARE  FOR  CHANNELS  11-15  IN  THE 
4.3  MICRON  BAND.  TO  USE  THIS  PROGRAM  FOR  ANY  SELECTED  CHANNELI S) 
CHANGE  THE  FOLLOWING  INITIALIZATIONS***** 

1. )  FILL  ARRAYS  "CNTR"  C "BWDTH"  WITH  DESIRED  CHANNEL  CENTER 

WAVE  NO.  AND  BANDW IDTH, RESPECT IVELY. 

2. )  MOLECULE  ARRAYS  I H2 0 ,CC2 , C3 , AN20, CO, CH4 , t02 ) SHOULD  BE  DIM- 

ENSIONED SO  THAT  THE  FIRST  ELEMENT  REFLECTS  THE  HIGHEST 
MOLECULE  NUMBER  DENSITY  FOR  ANY  CHANNEL  BEING  CONSIOEREO. 
THE  SECOND  ELEMENT  REMAINS  AT  4. 

3. )  ARRAY  " MOL 3 NO"  SHOULO  REFLECT  THE  SAME  UPPER  BOUND  AS  DE- 

SCRIBED IN  NO.  2 ABOVE. 

4. )  THE  DATA  IN  EACH  MOLECULE  ARRAY  BY  COLUMNS  IS***** 

1  WAVE  NUMBER 

2  LINE  INTENSITY 

3  HALF  WIDTH 

4  ENERGY 

FOR  EACH  CHANNELI S)  TO  BE  READ,  ONLY  THOSE  MOLECULES  PRESENT 
WITHIN  THE  BANDWIDTH  APE  WRITTEN  INTO  THE  MOLUCULE  ARRAY! S ) • 

THE  MAIN  PROGRAM  WILL  NOT  PRINT  A MOLECULE  HEADING  IF  THERE  ARE  NO 
MOLECULES  FOUND  ON  THE  TAPE  BETWEEN  THE  SPECIFIED  LIMITS. 

BEFORE  EACH  CHANNEL  OUTPUT  THE  CHANNEL  NO., CENTER  WAVF  NO.,  ANO 
THE  CHANNEL  BANOWIOTH  ARE  DISPLAYED.  AF  TER  EACH  CHANNEL  OUTPUT, 
THE  ENERGY  CONSTANT  AND  INTENSITY  AVERAGE  FOR  THE  CHANNEL  ARE  ALSO 
DI SPLAYED. 


INI TI ALIZATI CNS 

DIMENSION  CNTPI5 ), ICNTI7) ,SAVE( 7),EAVEI7» 

DIMENSION  ISUMI7)  ,BW0THI5) , MCLBND IT  I , T EMOL !4 1 

GOIMON  H20I  85,  4), C02I  1415,41,031  1,4), AN20I  800 ,4)  ,COI  20, 4)  , 

1 CH4I 1 ,4) ,0211,4) 

'DIMENSION  IFILE15) 

DATA  IFILE/81,82,83, 84,85/ 

DATA  MOLBND/85, 1415  ,1  ,800,20,1,1/ 

DATA  CNTR/ 2 190. 0,22 12.0, 2242.0,2275. 0,2357. 0/ 

DATA  I SUM/7  *0/ 

DATA  BWDTH/  20.6,22.5,21.6,35.2,23.0/ 

DEFINE  FILE  81(1000,  16,L,!A1) 

DEFINE  FILE  82  1 1 075 , 1 6 , L , I 82  1 
DEFINE  FILE  83 ( l 000, 16, L , IC31 
DEFINE  FILE  84  1 1 6 DO, 1 6, L , ID4  ) 

OEFINE  FILE  85(550, 16,L,IE5) 

JR-0 

SETTING  UP  CHANNEL  COUNTER  AND  CHANNEL  LIMITS 

DO  10  11-1,5 
JR= JR ♦ l 

»,  BESTfAVAILABLE  COPY 


1 


FOR  TRAN  IV  G LEVEL  21 


MAIN 


DATE  - 77108  11/31 


0018 

0019 

0020 
0021 
0022 


0023 

302* 

0025 


0026 

0027 

0028 

0029 

0030 

0031 

0032 
C033 
003* 
0335 
00  36 
0037 
0338 
0039 
00*0 
00*1 
00*2 
30*3 
00** 
00*5 
03*6 
00*7 
00*3 


JJMFILEI  ID 
ICHAN*  13  *■  II 

FHI  - CNTR(ll)  * I8WDTHI  ID/2.0) 

FLO*  FHI  - BWOTHt III 

WRITE  (6 ,100)  ICHAN,CNTP(  ID,  8WCTHI  ID 
C 

C CALLING  CHANNEL  REAO  SUBPROGRAM 
C 

CALL  CHREAOIFHI , F LO, MOLBND.JR , ICNT , SAVE ,E  AVE , I TOTAL) 
DO  111  L*l,7 

111  I SUM!  L)  = ISUMIL)  ♦ ICNTIU 
C 

C OUTPUT 
C 

ISTOP  = ICNT«  1) 

IF  I 1ST  OP . EQ .0 ) GO  TO  31 
WRI TEI6.ll) 

WRITEI6, 12) 

DO  200  120=1, ISTOP 

WRITEI6,!*)  I20,(H20II20,N) ,N«1,*) 

233  CONTINUE 

31  CONTINUE 
ISTOP=ICNT( 2) 

IF  I ISTOP. EQ.O)  GO  TO  22 
WRI TEI6.15) 

WR  ITEI6,  12) 

00  300  130=1, ISTOP 

WRITE  (6,1*)  I30,IC02(I30,N)  , N= l , * ) 

303  CONTINUE 

32  CONTINUE 

I S TOP= I C NT! 3) 

IF (ISTOP. EQ.O)  GO  TO  33 
WRI TE (6, 16) 

WRITE16.12) 

00  *00  1*0=1, ISTOP 

WR I TE ( 6 , 1 *)  1*0, (03(1*0, N) ,N=1 ,*) 

*33  CONTINUE 


00*9 
00  50 

0051 

0052 

0053 
005* 

0055 

0056 

0057 
. 00  58 

0059 
0C60 
0061 
0062 
006  3 
006* 

0065 

0066 

0067 

0068 
0069 


33  CONTINUE 

I STOP= ICNT ( *) 

. I-f  ((STOP.  EQ.O)  GO  TO  3* 

WRI  TE  ( 6 , l 7) 

”WRITF(6,12) 

DO  500  150=1, ISTOP 

WR I TE ( 6 , 1 *)  I 50, ( AN20 ( I 50 , N)  ,N*1,A) 

500  CONTINUE 
3*  CONTINUE 

ISTOP=!CNT( 5) 

IFI ISTOP. EO.O)  GO  TO  35 
WRl.TE  ( 6 ,*5) 

WR  I*TE(  6,  12) 

DO  600  160=1, ISTOP 

WR I TE (6,1*)  160, (C0(I60,N) ,N=1,4) 

600  CONTINJE 
35  CONTINUE 

ISTOP* ICNT( 6) 

IF ( ISTOP. FQ.O)  GO  TO  36 
WRI TE I 6 , *6) 

BEST  "AVAILABLE  COPY 
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MAIN 


DATE  » 77108 


11/S 


0070 
00  71 

0072 

0073 
00  7* 
00  75 

0076 

0077 
00  78 

0079 

0080 
0081 
0C82 
0083 
0)8* 

0085 

0086 


0087 

0088 
0089 
C090 

0091 

0092 

0093 
009* 
0)95 

0096 

0097 

0098 

0099 

0100 
0101 
0102 

0103 

010* 

0105 

0106 

0107 

0108 
0109 
01  10 
01 11 
0112 
0113 
Oil* 

0115 

0116 


C 

c 

c 


Of)  700  170=  1,1  STOP 

WRITE (6, 1*)  170, (CH*( 170 , N ) , N= 1 , * ) 

700  CONTINUE 

36  CONTINUE 

I STOP=I CNT (7 ) 

IFI  ISTOP.EQ.OI  GO  TO  37 
WRITEI6, *7 1 
WRI TEI6 ,12) 

DO  800  180= 1,1  STOP 

WRITE (6,1*)  180, !02(l80,N)»N=l*A) 

800  CONTINUE 

37  CONTINUE 
WRI TE(6 ,25) 

DO  65  19=1,7 

65  WR IT E (6,66)  19 , E AVE ( 19 ) ,S AVE ( 19 ) 

WRI TE (6, 112)  ( I SUM(L) ,L=l ,7) 

I STOP=0 

OUTPUT  PERMANENT  FILE  OATA  FOP  VERIFICATION 


WRITEI6.12) 

00  900  K= l , 7 

IFI  ICNT(K) .EQ.O)  GO  TO  900 

1 S TART  = 1 ST  OP  1 

I STOP= I STOP  «■  ICNT(K) 

DO  910  L=  1ST  ART,  ISTOP 
REAO(JJ'L)  (TEMCLINI ,N=l ,*l 
KK=L+ 1 
FIND(JJ'KK) 

WR  ( TE  ( 6 , 1 * I L,  ITEMOHM)  ,M=1,*) 

910  CONTINUE 
900  CONTINUE 

WRITE  (6, *2)  ITOTAL 
10  CONTINUE 

*2  FORMAT (IX,  • ITOTAL  IS', 110) 

112  FORMAT(///,25X,MCOUNT  TOTALS  FOR  ALL  CHANNELS  ARE  • ,/ , 23X  , 7 I 15,  i 


1 ) 

11  F0RMAT(1H1 ,25X,25HWATER  VAPOR  MOLECULE  DATA,/, 26X, 251 IH*)) 

12  FORMAT! 2X, 10HRECCR0  NO. , 7X , 8HWA VE  NO. ,3X ,1*HLI NE  I NT ENS I TV , *X, 
1V0HHALF  WIDTF,  7X , 6HENER  GY , / , 2X,  69(  IH- ) ) 

. 1*  FORMAT!  5 X, 15, 8X, FI 0.3 ,5X , El 0 .3 ,7X , F5  .3, 8X , F10 .3 ) 

25' FORMAT!///, 6X, 29HENERG Y AND  INTENSITY  AVER  AGES , / ,6X, 29 ( IH») , // ,! 

17HM0L  NO. ,5X,10HENESGY  A VE, 5X , 13HI NT ENS IT Y AVE, /, 5X, *01  IH? ) ) 

15  FORMAT!  1H 1 , 23X  ,2  BMC  AR  BON  OICXIDE  "OLECULE  .DATA,  / , 2 *X  ,28 1 IH* ) I 
IS  FORMAT! IH1.23X, 19H0Z0NE  MOLECULE  DAT A , /, 29X , 19 ( IH* ) ) 

17  FORMAT! IH1 ,2*X«c7HNITPCUS  OXIDE  MOLECULE  DATA, / , 25X, 27(  1H* ) ) 

*5  FORMAT!  IHI,22X,29HCAR80N  MONOXIDE  MOLECULE  D AT A ,/ , 23X , 29  C1H* ) ) 
*6  FORMAT ( IH1 , 2 7X , 2 1 FM ETHANE  MOLECULE  DATA ,/, 28X, 211 IH*) I 
47  FORMAT! IH1 ,28X  ,20HOXYGEN  MOLECULE  DAT  A,/ , 29X , 20 ( IH* ) ) 

66  FORMAT! 8X, I l , 8 X, F 10. 3 , 8X , E 10. 3) 

100  FORMAT!///, 5X ,1 AHCHANNEL  NO.  IS, I*,  3X* 25HFOR  CENTER  WAVE  NUMBER 
1,3X,F6. l,2X, 21HCHANNEL  BANDWIDTH  1S>,F6.1) 

STOP 

END 
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0031 


0002 

0003 

000* 

0005 

0006 

0007 

0008 

0009 

0010 
0011 


C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  CHREAD 
MICHAEL  T.  POTTER 

SUBROUTINE  "CHREAD"  WILL  EXTRACT  OATA  FROM  THE  "MCCLATCHEY"  TAPE 
BETWEEN  ANY  TWO  WAVENUMBER  LIMITS,  AND  BUILD  PERMANENT  OATA  FILES. 
OATA  IS  ALSO  bEGRAGATED  BY  PCLECULE  TYPE.  MOLECULE  COOES  ARE 
NUMBERED  FROM  1-7  AND  ARE  GIVEN  BELOW***** 

1.  WATER  VAPCR 

2.  CARSON  DIOXIOE 

3.  OZONE 

4.  NITRCUS  OXIDE 

5.  CARBON  MONOXIDE 

6.  METHANE 

7.  OXYGEN 

PARAMETERS  PASSED  TO  "CHREAD"  ARE  THE  HI  t LO  WAVE  NO  LIMITS  OF 
THE  DESIRED  BAND  AND  THE  ARRAY  CONTAINING  THE  UPPER  BOUNO  FOR 
MOLECULE  DENSITIES.  CHREAO  WILL  PLACE  THE  FOLLOWING 
DATA  IN  THE  FOLLOWING  ARRAYS***** 

1.  ICNT—  70  VECTOR  WITH  TOTAL  MOLECULE  NUMBER  DENSITY  FOR 
EACH  OF  THE  7 MOLECULES. 

2.  SAVE.EAVE—  7D  VECTORS  CONTAINING  THE  AVERAGE  INTENSITY  ANO 
ENERGY  FOR  EACH  MOLECULE. 

3. H20,C02,. ...  — MOLECULE  ARRAYS  WHICH  MUST  BE  IN  COMMON  WITH 

THE  MAIN  PROGRAM.  DATA  CONTAINED  IN  THE  * COLUMNS  ARE***** 

1  WAVE  NUMBER 

2  L INE  INTENSI  TY 

3  HALF  WIDTH 

* ENERGY 

THE  FOLLOWING  ERROR  FORMAT  WILL  BE  PRINTED  BY  "CHREAD"  IN  THE 
GIVEN  CIRCUMSTANCES***** 

1.  FORMAT  97-  DATA  FILE  ON  TAPE  DOES  NOT  CONTAIN  LOW  BOUNO  OF 

BANOWIDTH. 

2.  FORMAT  96-  THERE  ARE  MCRE  RECOROS  BETWEEN  LIMITS  THAN  THE 

MOLECULE  ARRAY! SI  WERE  DIMENSIONED  FOR  . 

AN  OPTION  EXISTS  TO  HAVE  THE  ROUTINE  PRINT  OUT  MOLECULE  OENSITY 
COUNTS  POR  EACH  TYPE  PRESENT.  SEE  COMMENT  CARD  AFTER  FORMAT  999 
ADDITIONALLY,  PERMANENT  FILES  81  - 85  ARE  PROOUCEO  AND  CONTAIN  ALL 
LINE  OATA  ARRANGED  BY  TYPE  IN  ASCENDING  WAVE  NO.  STAAT  ANO  STOP 
. POINTS  FOR  DATA  TYPES  CAN  BE  FOUND  BY  USING  ARRAY  ICNT.  EACH 
RECPRO  CONTAINS  FOUR  ELEMENTS  CORRESPONDING  TO  THE  COLUMNS  IN 
' THE  MOLECULE  ARRAYS. 

SUBROUTINE  CHRFADIFHI , FLO, MOLBND, JOE , ICNT  ,.SAVE,EAVE« I TOTAL  I 

CLEARING  MOLECULE  ARRAYS  ANO  INITIALIZING 

DIMENSION  WAV! 40)  , STRGH I *0)  ,WDTHI40) , ENGY (40 ) , ROT  I 0140 , 5 ) , V I BI40 > 
DI  MENS  ION  10 AT  E < 40 1,  ISO ( 40  I , MOL ( 40 J , GBGI 231 , ICNT! 7) 

01  PENSION  EAVE (71 ,SAVE(7) , KC LBNDI7  ) 

DIMENSION  IF ILEI 5) 

COMMON  H2 0 1 8 5 , 4 ) , C02 1 1 41 5 , 4 ) , 031 1, 4) , AN20 1 800 , 4 ) ,C0( 20, 4) , 

1 CH4I  1,  4)  ,021  1,4) 

OATA  MH20, MC02, M03*  M AN20, MC0,MCH4,M0  2/  85,1415,1,  800,20,1*1/ 

DATA  IFILE/81,82,83,84,85/ 

DATA  EAMAX, EAM  IN/O  .0, 1.0E10/ 

DO  87  MM*l ,7 
l CNT ( MM ) = 0 
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0012 
0013 
C014 
00  15 
0016 
OCl  7 
0018 
CCl  9 
0020 
0021 
0022 
00  ?3 

0024 

0025 

0026 
00  2 7 
0028 
0029 


0030 

0031 

0032 
00  33 
00  34 

0035 

0036 

0037 

0038 
00  39 
004  0 


0041 


0042 
00  43 

0044 

0045 

0046 

0047 

0048 


0049 
00  50 
0051 
005? 
00  53 

0054 

0055 
00  56 

0057 

0058 


C 

C 

C 


C 

C 

C 


C 

c 

c 


SAVE!  MM ) = 0.0 
E AV  E I MM ) = 0 • 0 
87  CONTINUE 

Dn  70  11=1,4 
00  71  K 4=  l , MH20 

71  H20( K4 , I I ) =0.0 
00  7 2 K 3= l » MCO  2 

72  C02(K3,II)=0.0 
00  73  KC  = l » M03 
03 (KC, I I ) = 0 .0 
CH4 ( KC , II ) =0.0 

73  02IKC, I I 1 = 0.0 

00  74  K 0=  1 , M AN20 

74  AN20I KO , I I 1 =0. 0 
00  75  K E=  1 , MCO 

75  COIKE, II)=0.0 
70  CONTINUE 

JO=IFILEl JOE) 

LOCATING  3EGINNING  WAVE  NUMBER  CN  APPROPRIATE  FILE 
NOF I L=l 0 

IFIFLO  .GT.0050C.0)  N0FIL=20 
IF (FLO  .GT. 01000.0)  N0FIL=30 
I F ( p LO  .GT. 02000.0)  NOF I L=40 
IFIFLO  .GT. 050)0.0)  NOF IL=  50 
IFIFLO  .GT. 07500.0)  N0FIL=60 
IFIFLO  .GT.10C00.0)  N0FIL=70 
05  RF  AO (NO  F I L , 9 1 , ENO=  99  ) INUM.FR 
IF( FR.tT.FLC)  GO  TO  05 
BACKSPACE  NOF  IL 
BACKSPACE  NCFIL 

READ  TAPE  FROM  LOWER  LIMIT  TO  HIGHER  ONE 

13  READ!  NOF  I L,  21  ,END  = 33  ) INUM,  (W  AV  I K ) , S T RC-H  (K  ) , W DT  H I K ) , ENGYIK  ), 

1 (ROT  I OIK, J) , J=l , 5) , VI  HI K) , IDATE(K)  ,1  SOI K)  , 

2 MOLIK)  , K=l,  40  ) , (GBG(M),,M=l,  23) 

00  20  12=1, INUM 

K-MOL  = MOL  I 12) 

IFIWAVI 12) .LT.FLC)  GO  TO  20 
• IF!  V.A  VI  I 2)  .GT.FHI  ) GO  TO  15 
ICNTIKMOL ) = ICNT I KMOL ) + 1 
J J=I CNT I KMOL) 

IF  (ICNTIKMOL)  .GT.  M0L8N0I KMOL) ) GO  TO  95 


DATA  SE3REG4TI CN 


IF ( KMOL. E 0.1)  GO  TO  101 
IFKM0L.E0.2)  GO  TO  102 
I F ( KMOL . E Q . 3 ) GO  TO  103 
IF< KMOL. EG. 4)  GG  TO  104 
IFIKMOL.EQ. 5)  GO  TO  105 
IF (KMOL. EO. 6)  GO  TO  106 
GO  TO  107 

101  H20 ( J J , l ) =WAV( 12) 

H?0(JJ,2)=STRGH( 12) 

H20( J J , 3 ) = WO THI12) 
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0097 

0098 

0099 

0100 
0101 
C102 
0103 
010* 

0105 

0106 
Cl  07 
0108 


0109 

0110 


' GO  TO  909 

99  W»ITP(6,97>  NOFIL.FR 
REWIND  NOFIL 
GO  TO  999 
33  REWIND  NOFIL 

IF1N0FIL.E0.70)  GO  TO  99 
NOF  IL  = NQF  IL  ♦ 10 
GO  TO  13 

95  HR  7 TE 1 6,96)  1 ICNT1 J)  , J = 1 ,7) ,WAV (12) 
REWIND  NOFIL 
999  CONTINUE 

WRITE! 6,  755)  1 I CN T( 66 ) , LL=l , 7) 

C 

C AVERAGING  ENERGY  AN9  INTENSITY 
C 

00  172  19=1,7 

IF( ICNTI 19) .EO.O)  GO  TO  55 
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0111 

0112 

0113 

0114 

0115 

0116 

0117 

0118 


0119 

0120 
0121 
0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 
01  40 
0141 


0143 

0144 

0145 

0146 


SAVE!  19  )=  SAVE  I 19)  /I  CNT I 19) 

EAVE ( 19)  = EAVE(  I9)/ICNTI 19) 

GO  TO  172 
55  SA  VE  I 1 9 ) = 0.  0 
EAVEI  19  ) = 0 .0 
172  CONTINUE 
K=  1 

ITOTAL  = 0 
C 

C BUILD  PERMANENT  DATA  FILES  FOR  CHANNEL 

C 

DO  500  150=1.7 

IF  (ICNTII50)  .EQ.O)  GO  TO  500 
KK= I CNT (150) 

DO  550  M=1,KK 

GO  TO  ( 561, 562 , 563,  564,  565,  566,  567)  , I 50 

561  WRITEIJO'K)  (H20  IM,N),N=1,4) 

GO  TO  550 

562  WRITEIJO'K)  (C02  IM,N),N=1,4) 

GO  TO  550 

563  WRITEIJO'K)  103  IM,N),N=l,4) 

GO  TO  550 

564  WRITEIJO'K)  ( A N20 ( M , N)  , N = l ,4) 

GO  TO  550 

565  WR I TE ( JO*  K)  ICC  IM,N),N=1,4) 

GO  TO  550 

566  WRITEIJO'K)  (CH4  (M,N),N=1,4) 

GO  TO  550 

567  WRITEIJO'K)  (02  IM,N),N=1,4) 

550  K = K + 1 

I TO  TAL=  ITOTAL  «■  ICNTU50) 

500  CONTINUE 

21  FORMAT! 1 10, 40  IF  10.3,610. 3,  F5 .3 , FI  0 .3, 5A6 , A5 , 13, I 4, l 3) , 23A 30) 

91  FORMAT!  T 10.F10.3) 

96  FORM  AT  I 1H1,23HCIMENS ION  SIZE  EXC EEDEO , / , IX , 7! I 6, 2X) , / , 1 X, 

1 20HLAST  WAVE  NUMBER  I S : ,2X , F10 .3 ) 

97  FORMAT!  1H2, 5X,9H  FILE  NO : , 2X , I 2 , 2X , 2 5HD0E S NOT  CONTAIN  WAVE  NO: 
12X.F10.3) 

755  FORMAT! IX,' ICOUNT  FOR  ABOVE  CHANNEL  NO,  IS • , / , 7 1 18 , 2X ) ) 

RETURN 
' END 
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XXSVSL IN 
XX 

//FOPT.SYSIN  01 


00  OSNAMF  = 'LOAD  SET, II  SP=(  POD, PASS)  .UNIT-SVSSQ, 

SPACE  = 180, ( 23).  I ) )), RLSE ) , OCB=3LKSIZE=8D 

t 

STEPNAME 


X 28000020 
*05019)2) 


IEF236I 

ALLOC 

. CDR  J03N  78 AO 

1 FORT 

1 E c2 T 7 I 

AA3 

allocated 

TO 

SYSPR  INT 

!?<=  2371 

431 

allocated 

TO 

SYSPUNCH 

\ CF237  I 

232 

ALLOCATE, 1 

TO 

SYSLIN 

! E c 2 3 7 1 

All 

ALLOCATED 

TC 

SYS  IN  - 

1 EE  1421 
!=F285  I 
I F c 2 3 51 


XXG.D 

xxsvsi.  r* 

xx 

xx 

XX  SV  SLP.iT 
XXSYSLIN 
XX FTO  5C )) l 


STEP  WAS  EXECUTED  - CONI  CODE  0000  - 

SYS  77  l )8 .T 165922.P  V0  1 1. J08N  7840. LOAOSE T PASjfb 

VOL  SE * NOS  = 'ITCPSA. 

J T E 3 1 1 1 STEP  FOPT  E XCC  l T I TIME  = 010.17  SEC. 

UTEP11I  STEP  F?RT  REGION  = 0123K  MAX  CORE  USEO  * 0102K 

EXEC  ?3K  = L0A.1ER,PAPM=(  WA  P , LET  , PR  I KT  ) ,CONO*  I A,  LT,  FORT) 

DO  9SVAWE-SYS  2 .FOR  TL I 3. DI SP=SHR 

10  DSN= DTEP.FQRTLI 3,0ISP=SHR  t FORTGCLO 

30  0 SN  = :JTEP  • C 4LC CM® . LOAD  , DI  SP-SHR  FORTGCLD 

DO  SYS  3U  T = A 

DO  DSN4M==*.FGRT.SVSLI N,OIS®=(  OLC.OELFTE ) 

OD  DON  AN  E = SYS  IN 

//3P.FT  16  F 1 1 1 DO  SYS  OUT -IA..SA10I 
X/CT06C001  OD  SYSOUT  = 4 

XXFT07C1DI  Dll  SVS1UT=3 

//GO. =T l OFOOl  ID  IIN|T=2400,OSN=CATA10,LA8FL=(02,BLP,, IN). 

//  DI S®=nLD,  V0L  = SC»  = T90339,DC3  = < RECF,M=F  .3LKSI ZE =39001 
//G0.F7 2 1F1 11  DD  UN  I f = ?4 1 1 , OS  N=0  AT  A 2 3,  L AB  EL  = I 15.BLP. . INI • 

//  DI  S®  = .1LD  ,V>  L --  JER  = TOD383,0C3  = (RECcv  = F,BLKS  IZE=3900  ) 

/ / GO  ,cT3DcDll  DO  UN l T= 2419, D SN=D4 TA 3 ), L ArtEL=C 08 »BLP » » I N)  • * 

//  DI ;p  = DL"  . /DL  = S?R  = T 11333. 0C9= I R ECF  M-  E t BLK S I ZE= 3000  I 
/ /GO.  c T4DC 1 1 l 1 . UNI  T=2  A 10.0  SN=0  A TA  AO  , L ABEL  = C 1 1 .8LP..IN). 

//  r»Isn=OLO.VOL=SFR=T00333,OC3=|RECFM  = F,3LKSIZE  = 3999)  * 

//G  3.  c T5  9FC01  O'*  UN IT =24 10, OS N =C 4T A50.LA6 EL =114, BLP,» IN). 

//Of  S®  = )LD«  VDL  = 3ES  = T99338, DC®  = I RFC e W=F  ,«LKSI  ZE  = 3 900) 

//GO.PTs  J Fill  Tj  UN  IT=24D1,0SN=0ATA6D,LA8EL= I 17.BLP..IN)  , 

//  II  S°=DLr> . VOL  = SER  = T0O3  3S  ,nC  8 = I R ECF  P = F , 3LKS  I ZE  = 3900  I 
//G0.PT7DF1DI  DO  UN  I T = 24 10. 0 SN =D A T A 70 , L AP EL  = < 20 , 3LP , , I N)  , 

//  01  SP=OL".7OL  = SER=T00:?33,0C9=  tREC^F,  3LK  S I Z E=  3900  I 
//GC.FT81PP DI  DO  ilSN=USE  R.  EFEO  294  7 . POTMOL . 0 1 S P= ( NEW , C AT  LG, DEL ET  E I . 

//  VCL=SER  = UTE°33,UN  IT=23X4, SP  ACE=( CYL , I 02,011 .RLSE ft 
/ /GO. F T R2FCQ1  00  DSN=USE 3 . EEF12S A7 . POT 833 , 0 IS P= (NEW , C ATL G, DEL ET E ) , 

//  V0L  = SER*'JTEP33,';';IT  = 2314,SPACE  = ICYL,(02,0l!  .RISE) 

//3D.tT33FG ii  oi  OSN= ilSFR  .EE ED  29 A7. POTCCC , 0 1 S®= (NE W, CATL G .DELE TE ) , 

//  V?c  * Sc  R = 'JTC  P3  3 . UN  J T=2314»SP AC E=ICYL, (02,01  ),RLSE) 

//G0.PT3AF1DI  OD  OSN  = USER .EEE020A7.o0T000.0l S°  = I NEW, CATLG .DELETE!  , 

//  VOL  = SP#  =IJTC  ®3  3,iJNIt  = 231A,SPACE=CCYL,  ( 12, 01). RLSE) 

//GO.cT85FOOI  OD  0SN=USER.EEE029A7.  P.DTEEE  , 1)1  SP  = I NEW  , CAT  LG , DELETt!  , 

//  VOL=SFR=jtepu,uN  !T=231A,  SPACE  = (CYL,(  01,01)  .RLSEI 
//GO. SYS!  N OD  • 

// 


42310120 

A9131D21 

49000035 

49100040 

56000020 

63000020 

709)0121 

77110020 

84)10)20 


IFF253I 

J03N78AO  GO  FT81F311  OIRECT  ACCESS-DUPLICATE  NAME 

I CF  23  5 1 

SYSL.F0RTLI9 

KEPT 

!=F2B5l 

VOL  *SER  NOS=  LNK2I8. 

IEF235  | 

UT  EP . FDR  TL 1 3 

KEPT 

I EF  2 351 

VOL  SER  NOS=  UTEP11. 

IEF285I 

UTEP.CALCOWP.LDAO 

KEPT 

1 EF285 I 

VOL  SE®  NOS  = UTEPll. 

IFF  28  51 

SYS  77108. TC6  592  2. R V001. J03N78A0.  LCADSET 

DELETED 

IEF285  I 

VOL  SER  NOS  = UTEP3A. 

IE  F 2851 

DATA  10 

KEPT 

!CF285I 

VOL  SER  NO S=  T00388. 

IFF285I 

OAT  A20 

KEPT 
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FOR  TRAN 


0001 

0002 

0003 

0004 


0005 

0006 

0007 
CC08 
00  09 
0010 
0011 
0012 

0013 

0014 

0015 
TO  1 6 
0017 
00  18 
0019 
00  20 
0021 
0022 


0023 

0024 

0025 

0026 
002  7 


0028 


0 029 
00  30 

0031 

0032 

0033 

0034 
00  35 


IV  G LEVEL  21 

■MAIN 

DATE  = 77108 

1 1/01/56 

C 

C 

ABSORPT ION 

COEFFICIENT  TEST  ROUTINE 

C 

C 

MICHAEL  T.  POTTER 

C 

C MAIN  ROUTINE  INITIALIZATIONS 
C 

REAL  *8  A BMCD ( 3 3, 7)/ 23 1*0. 00000000/, S TDK!  33,  71/231*0.00000000/, 

1 PEPDIF I 7) /7*0.  00000000/ 

REAL *3  A3DIFI 71/7*0.00000000/ 

REAL  *3  P(33) , PR (33 ) ,T (33) 

R F A L * 8 SPEC  1 ( 1600) 71600*0. COCOOOOO  / , SP tC 2 < 1 600  1/16  00*0.00000000/ , 

1 SPEC3< 16  TO)/ 16 0 0*0.  TO  TOO  TOO/, SPFC4<  1600) / 1600*0.  00000000/, 

2 SA  VI T(  7 , 9,71/441*0.00000000/ 

REAL  * 8 AA,ALFA,AMARG,63,3LTZ  ,C  N VRT  1 , CN  VR  T.2  , ENG  , ONE  , P I E , P I V , 
l °LK, POV, PO, SPCL , TEST ,T0,WA,WA2 
REAL  *8  ENSA ( 7 1 /2. 936C8275E-20, 5. 73070102 E-20, 0.00000000 , 
l 2 .28626414E-20.  1 . 7 55  3 19  3 2E- 2 0, 0. CO 000  TOO , 0.  00000000/ 

DIMFNSICN  HGT ( 33  > , I CNT ( 7 1 , CNST G 1 7 1 , MOLNAM ( 7 ) 

0 I ME  NS  I CN  Nc  W l ( 7)  ,NEW2(7)  , Ft (100)  ,FA(100)  ,FW1100) ,FS(100) 

DATA  LEVEL , ITOT AL/33, 593/ 

DATA  CNSTG/ 1.5, 1 .0,1 .5 ,1 .0 ,1 .0, 1 .5, 1 .0/ 

DATA  N=Wl,NEW2/7*0, 7*0/ 

DAT  A 8L  TZ  , PLK,  S» OL/  1 . 3 306 2 200E-23, 6 . 6 26 l 9600E-34, 2 . 99792 50 02  *-10/ 
DATA  WA  , TO, PCI, PIE  ,AA,B3/l.  18400000,2.96000000  202,  l .003  00000, 

1 3 .141592  65, .65TTT0TO, . 70TCC0T0/ 

DATA  ONE , A MARG  , VCNT , NCNT / 1 .C3000000 , 0 . 0 52 50000 , 7 , 9/ 

DATA  IC  NT / 34,195,0,743,16,0,0/ 

DEFINE  FILE  81  ( 2 50 3,  16, L , I A ) 

WAV-2  190. 0 

CNVRTl=SPOL*PLK 

CNV RT2  = ONE/ ( RLT  Z*TO ) 

PI V=ONE/Plc 
POV=ONF/PO 
WA2  = ONE  -WA 
C 

C INPUT  ATMOSPHERE  PROFILE 
C 

DO  10  I 1 = 1 .LEVEL 

READ!  5,01)  HGT(  T1),PR(I1)  , T ( I 1 ) 

C CONVERT  PRESSURE  FROM  MILLIBARS  TO  ATMOSPHERES 
P( I 1 ) =PR ( I l ) *9. 872-04 
10  CONTINUE 

RE  AO ( 5 , 32 ) I MCLNAM(k) ,K  = 1 ,7) 

C 

C COMPUTATION  OF  STANDARD  LORENTZ  LINE  ABSORPTION  COEFFICIENT 
C 

CALL  AB  STD  I I CNT , WA V , L E VEL , T , P ,C NS TG , I T OT AL ,ST DK) 

C 

C OUTPUJ  CF  COEFFICIENTS 
C 

WR I T E (6 , 03 ) 

WRITE  I 6, 04) 

WRITE(6,05)  (M0LNAM(K),K=1,7) 

DO  30  I 3=1 .LEVEL 

WR  ITEI  6,  06)  HGTI  I 3)  .(  STDM  I 3 ,N)  »N=  1 , 7) 

30  CONTINUE 
WR  ITEI  6, 19) 


BEST  AVAILABLE  COPY 


r 


F CP  T 3 AM  IV  G LEVEL  21 


C 

C 

c 

c 


COMPUTATION  OF  MODIFIED  LORENTZ  LINE  ABSORPTION  COEFFICIENT 


0036 

IHALT=0 

0037 

Nl  =0 

0038 

N?  =0 

00  39 

DO  90  19=1,7 

0040 

IF11CNTII9)  .FQ.  0)  GO  TO  90 

0041 

I G0=1  + [HALT 

00  42 

IHAL  T=  I HAL  T + ICNTI I 5) 

3043 

DO  100  113=  IGO,  I HAL  T 

00  44 

READ! 81  * 1 101  WNUM.ST ,ALFAC,EY 

0045 

JK=  1 10  ♦ 1 

0046 

FINC181 ' JK) 

0047 

TEST=ONE*(  WAV-WNUM) 

0043 

IF ( DABS ( T EST ) .LT . AMARG)  GO  TO  101 

00  49 

ALFA=OME/ALFAO 

00  50 

ENG=CNVRT1*EY 

0051 

Nl = Nl  + l 

0052 

SPEC11N1) = TE  ST* ALF A 

3053 

SPEC2IN1)  = (ENG •ENG At  I 9) ) *CNVR  T2 

0054 

SPEC3INI)  = ALF A*ST*DFXP ( ( ( ENG*WA2 ) + ( ENGA ( 15  ) *W A )) *CNVRT2) 

0055 

SP  EC  4(  N l ) = ONE/(  SRECUNl)  *SPEC  1 INI  1 + AA) 

0056 

NFWl  ( 19 ) = NEW  1 ( 19  ) + l 

0057 

GO  TO  100 

0058 

101 

N 2=  N2  + l 

0059 

NEW2 l 19 ) = NEW? (19)  + 1 

00  60 

FW( N2)=WNUM 

3061 

FS(N2)=  ST 

0062 

FA( N?)  =A  LF AO 

0063 

FE ( N 2 ) = EY 

0064 

100 

CONT  INUE 

0065 

90 

CONTI NUE 

0066 

I OU  I T=0 

0057 

DO  200  120=1,7 

0063 

IF(MEW1( 120) .EQ.O)  GO  TO  200 

0069 

ISTART=1  «■  IQUIT 

0070 

I OU I T=  IQUIT  + NEWl (120) 

0071 

DO  300  130=1, NCNT 

007? 

'NPO=J30  - l 

OC73 

DO  400  140=1, MCNT 

0074 

M=  140  - 1 

0075 

IF ( ( 40  .EQ.  MCNT ) GO  to  50 l 

00  7 6 

DO  500  I 50  = 1 START, I OU I T 

0077 

SAVIT( 140, 13 3, 120)=  SA VI Tt I 40 , I 30, I 20)  + ( SPEQ4 ( I 5 0 ) ** l 40 ) * 

l 

(SPECS II50)**NPC)*SPEC3(I50) 

00  78 

500 

CONTINUE 

0079 

GO  TO  430 

0030 

501 

DO  50?  152=ISTAPT, IQUIT 

0081 

SAVITI  I 40,1 3C, I 20)  = SAVI T( 140,130,120)  ♦ (SPEC4  (!52)**M)* 

1 

(SPEC2 (152 !»*NPO )*(SPEC3(  152))/(SPEC1(I 52)* SPEC  1( l 52) ) 

0082 

502 

CONTINUE 

0083 

400 

CONT  INUE 

0084 

300 

CONTINUE 

0035 

200 

CONTINUE 

0086 

CALL  ABSM00(NEWl,NEW2,T0,P0V,T,P,LE VEL , A A , 8B , MC NT , NC NT ,E NG A 

1,PIV,FW,FS,FA,FE, S A V l T , CNS TG , W A , WAV , ABMGD ) 


BEST  _ AVAILABLE  COPY 


125 


FORTRAN 

IV  G LEVEL 

21  MAIN  DATE  = 77108 

0087 

WRI TE (6 ,03) 

* 

0089 

WRITE! 6, 07) 

0099 

WRITE (6, 05)  (M CL NAM (K ),K=1, 7) 

0090 

00  40  14=1, LEVEL 

0091 

WRITE! 6, 06)  HGT ( I4),( A3M00!  I4,M)  ,M  = 1,7) 

0092 

40 

CONTINUE 

0093 

r 

WRI TE! 6,19) 

L 

c 

r 

CALCULATION  AND  OUTPUT  OF  PERCENT  RELATIVE  ERROR 

0094 

WR  ITE(6, 39) 

0095 

WRI TE (6 .1 l ) ( MCLNAM(M) ,N  = l ,7) 

0096 

DO  50  I 5=1, LEVEL 

3097 

DO  60  16=1,7 

0093 

I F< STOK! I 5, l 6)  .EQ.  0.30000000)  GO  TO  60 

0099 

R ER  D I F ( 16)  = 10  3.3»DA3S( ( STOK! 15,1 6)-A3M0D! I5,I6))/STDK(I5 

01  00 

60 

CONTINUE 

0101 

W° ITE(6,06)  HGT! I 5) , (PEROIF(L) ,L=1 ,7) 

01  02 

53 

CONT  INU6 

L 

c 

CALCULATE  ABSOLUTE  ERROR  IN  INVERSE  A TM— CMS 

11/01/56 


0103 

0104 

0105 

0106 
0107 
01  03 
01  00 
0110 
0111 
0112 

0113 

0114 

0115 

0116 

0117 

0118 

0119 

0120 
0121 


WRI TP(6 ,18) 

W«UT5<  6,  1 l ) { MOL  NAM  ( MMM)  ,MMM=l ,7) 

00  51  151=1, LEVEL 
DO  61  161=1,7 

6 1 ABOIFI  161 )=2.69E+  19*0 ABS(  STOK(  I 5 1 , I 6 1 ) -4  3 MOO ( 151,161)) 

51  WRITE  (6,06)  HGT< I5l>,  CAROIF(  I ),  1=1,7) 

01  FORMAT! E 5. I , 1X,E9.3,F 5. 1) 

02  FORMAT (7 A4) 

03  FORMAT! IH1 , 52X ,26HABS0RPT ICN  COEFFICIENT  FOR) 

04  FORMAT!  IH  , 5 l X , 2 7HSTANDAR0  LORENTZ  LINE  SHAPE) 

05  FORMAT !1H  , 5 1 X , 28 ( 1 H* ) ,// , 3X , 7HA LT ( KM ) , IX, 7 ( 5X , 44 , 6X ) , / , 1 X , 1 1 5 
1! 1H-) ) 

06  FORMAT ( IH  , 2X  , F 5 . 1, 3X , 7!  IX , F 1 3. 7, IX) ) 

07  FORMAT ( 1 H ,5 IX , 27 HMOD I F I ED  LORENTZ  LINE  SHAPE) 

09  FORMAT! IH l, 54X, 22HPERCENT  RELATIVE  ERROR ,2 ( / ,55X ,22 < IH*) ) ) 

11  FORMAT ( IN  ,2X,7FALT(KM>, IX , 7 ( 5X, A4 , 6 X ) , / , IX,  115!  1H-) ) 

1 8 FORMAT! IH1 ,43X .33HA9SCLUTE  ERROR  IN  INVERSE  AT M-CMS, /, 2C 49X, 331 IH* 
1),/)) 

19  FORMAT!//, IE  ,441-UNITS  ARE  IN  SQUARE  CENTIMETERS  PER  MOLECULE) 

' STOP 

ENO 


126 


BESOVAIWBLE  COPY 


FORT0  AN  IV  C LEVFC  21 
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0001 
00  02 

0003 

0004 

0005 

0006 
0007 
0C08 


00C9 

0010 

0011 

0012 

0013 

0014 


0015 

0016 
0017 
CC18 
0C19 
0020 
C021 
00  22 

0023 

0024 
C02  5 
0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 
00  36 
0037 


C SIJ  BRC'JT  IN  c ABSTD 

C MICHAEL  T.  PCTTER 

C THIS  SUBROUTINE  USES  THE  STANDAPO  LORENTZ  LINE  SHAPE  ANO 

C COMPUTES  ABSORPTION  COEFFICIENTS  L 1 N E-SY-L IN E . 

C COEFFICIENTS  ARE  PRODUCED  FOR  EACH  GAS  AT  EVERY  LEVEL.  THEY  ARE 

C PASSED  BACK  TO  THE  MAIN  SOUTINE  IN  ARRAY  STDK. 

C 

SUBS  OUT  I Nc  A3  STD ( I C N T , C TRFR 3 , LE V EL , T FMP , °RES S , CNSTG , I TOT AL ,STDK1 
IMPLICIT  RRAL*8I A-H.O-Z ) 

REAL  CNSTG.CTRFPC , ^AVO.STRO, ALFAC.ENGO 
DIMENSION  AL  FAD ( 1600)  ,ENG0(  1600)  ,WAV0(  16  00)  .STROU6001 
DIMENSION  ICNT (7 ) ,STDKI 33, 7) ,QRTC( 7) , TEMPI 331.PRESSI 33J .CNSTGI  71 
DATA  8LTZ * aLK , S POL/  1 . 38 062200 E-23 , 6 . 6261 9600 E-34, 2 . 99792 503E  + 10/ 
DATA  PIE,T(0,PO/3.  1415926  5,  2.  96  C SO  300E  02 , 1 . 0 0 009000/ 

PI EF=l .0/ PIE 
C 

C ATMOSPHERIC  LEVEL  LOOP 

00  5 15=1  , I TCT AL 

READ! 81 M 5)  WAVOI I 5)  , STH 0 ( I 5)  , ALFAOI 15)  ,ENGO( 15 ) 

JA= 15+1 
FI  ND(  81»  JA) 

5 CONTINUE 

00  10  1 1 = 1, LEV  EL 
C 

C ONCE  PER  LEVEL  CALCULATIONS 
T1=TP/TFiMP(  II  ) 

R 1 = PRE iS( III /°0 
RT0=P1*0S0RTIT l) 

DFNCM=( 1 .0/ I 8LT  Z*TO) ) 

DO  15  I 5=1,  7 

15  STDK! II , 15 >=0.00000090 
TP  =TF  MP ( I 1 ) 

CALL  RARTI TP ,QRTO) 

1ST  OP=0 

C MOLECULE  GAS  LOOP 
DO  20  12=1,7 

1 f ( I CNT I 12) . EQ .0 ) GC  TO  20 
I STAR  T= l + I STOP 

l S T 0 P = ISTOP  + ICNT I I 2) 

TPVJR  =Tl  **  CNSTG  I 12  ) 

\?"l3  = 0STn(  12) 

DO  300  13=  ISTART,  ISTOP 

STDK  I I 1 ,1  2I=STDK(  11  , 121  + TEMPKt  WAVOI  1 3 ) , S T RO  I I 3 1 , ALFAOI  13  ),  ENOCH  l 
13) , OFNOM , TP WR , T i »R TO, C TRFR Q) 

300  CONTINUE 

STOKI I 1, I21=STDK( 1 1 ,12) *VI 3*  PI  EE 
20  CONTINUE 
10  CONTINUE 
RETURN 
ENQ. 


BEST  AVA1UBIE  COPY 


. 


FOR  TRAN 


OOOl 


0002 

0003 

0)34 

0005 

OC  06 
003? 


0008 

CC09 
0010 
3 311 
0012 
3013 
C01  4 


0015 

0016 
0017 
00  18 

00 10 
0020 
C021 
0022 

0023 

0024 
3025 
0026 
0027 
00  28 

0029 

0030 

0031 
00  32 

0033 

0034 


IV  G LEVEL  21 


MAI  N 


DATE  = 77108  U/Ol/56 


C 

C 

C 

C 

C 

C 

c 

c 


c 

c 

c 


c 

c 

c 

c 


c 

c 

c 


c 


c 


SUBROUTINE  ABSMOO 
MICHAFL  T.  POTTER 

THIS  SUBROUTINE  CCut>UTES  AN  ABSOFBTION  COEFFICIENT  USING  A 
M30IFIE0  LOPENTZ  LINE  SHAPE.  COEFFICIENTS  ARE  PRODUCED  FOR  EACH 
MOLECULE  TYPE  EVERY  ATMOSPHERIC  LEVEL. 


SUBROUTINE  4R S MOD ( NE W l , NE W2 , TO ,POV,T  ,P  , LE VE L , A A , 68 , MCNT , NCNT , ENGA, 
ICNVRT2,  PIV.  FW,  FS, FA, FE, SAV  IT  , CN  STG  , Vs  A , WA  V,  AdMOD  I 


INITIALIZATIONS 


IMPLICIT  RFAL*8< A-H.O-Z) 

REAL  WVUM,  ST,EY,  ALcAO,CNSTG,FV»,FS,FE  , FA,  WAV,  FAC 

RE AL*3  CHEBYI9 )/ 1 .00000000, .99997000, .49999100,.  166  76500, 

1 . 41683 400E-1 , . 30324750 ,. 00137756, .2 2577900E-3, .27 374800 E-4/ 
DIMENSION  NEwl 1 7) ,NEW2( 7) , SAVI T( 7,  9 , 7) , FE I 1 301  , F S 1 1 3 3 J ,F A 1 1 )3» , 
1 FWI 100 ) , CNSTG ( 7 ) ,T 133) , P (33), ABMODI 33, 7 ) , OR  TO ( 7 ) , F AC  I 1 1 ) 
DIMENSION  ENGA I 7) 

DATA  FAC/ 1 .3, 1.0, 2.0,6. 3, 24. 0, 120.0, 720. 0, 5340. 0 ,403 20.0 , 

1 362330.0,3623300.0/ 


ATMOSPHERIC  LEVEL  LOOP 


DO  200  1 20  = 1 , LEVEL 
CALCULATE  PARTITION  FUNCTION  VALUES 
CALL  PART( TI 120) , QRTO) 

rcNST=To/r<  reel 

PCNST=P( 120 )*POV 
WA  ta  = TCNST*CNVRT2 
X=PCNST*OSQRT( TCNST) 

GM=l A4-(X*X  ) ) 

ABSORPT ION  GAS  LOOP 


502 

503 
400 


DO  250  125=1,7 

IF(NCW1II25)  . EQ  .0 ) GO  TO  253 
SAVEX=O.OOOCOOOO 

UAT^P !V*QRTO( I 25) *( TCNST**CNSTG I 12  5) 1 *DEX P ( -ENG A I I 25 ) * WAT A ) 

N COUNTER  LOOP 

DO  300  130=1, NCNT 

NPO= 133-1 

SAVGEf =0. 09000000 

E XPON=  I < WA-TCNSTJ^^NPO^CHEBYI  130) 

M COUNTER  LOOP 
DO  400  I 40  = 1 , MCNT 


M= I 40-1 
GEE= ( GM**M ) *X 

I c I'M  4 0 .EC. MCNT)  GO  TO  502 
S A V 1=  SA  V 1 T t 140,  130,  125) 

GO  TO  503 

GEE  = 3B*GFR  / ( X*  X*  BB) 

S AV I =S AV 1 T I 143,  130,  125) 
CONTINUE 

S AVGEE=  S AVGEE  *■  SAVI*GEE 
CONTINUE 

SA  VE  X = SA  VE  X*  EXPON*SAVGEE 


128 


B®  WMVABll  COPT 


FORTRAN  IV  G LEVEL  21 


4BSM0D 


DATE  = 77108 


11/01/56 


300  CONTINUE 

ABMODI 123, 125)=  WAT«SAVEX 
250  CONTINUE 
255  1 S T0P  = 0 

03  45  1=1,7 

H0L0=0. COCOOOOO 

IF  INEW2I I ) .EC.  0)  GO  TO  45 

IBEGIN  = l + ISTOP 

I STOP  = ISTOP  ► NEW2 ( l ) 

TPW=TCNST**CNSTGU  ) 

03  35  I3=13EGIN* ISTOP 
WNUM=F  wi  I 3) 

st=fsi m 

ALF  A0=c  A I 13) 

EY  = FF(  r 3) 

STO=TENPK  (WNUM.ST,  Al.FAO,  EY,CNVRT2,TPW,  TCNST.X,  WAV) 
H3L0=H0L3+$TD*CRTCI I )*PIV 
35  CONTINUE 

ABMODI 120,  I )=ABMOQ(  120,  I ) ♦ HOLD 
45  CONTINUE 
203  CONTINUE 
RETURN 
ENO 


BESOVAILABLE  COPY 


POFTOAN 


0031 

0002 

0033 

0004 

0005 
00  36 
0007 

0009 
OOC9 

0010 
0011 
0012 


IV  G LEVEL  21 


MAIN 


DATE  = 77108  11/01/56 


FUNCTION  SUBPROGRAM  FCR  THE 
COMPUTATION  pF  A8S0TPTICN  COEFFICIENTS 
USING  STANC4RD  LORENTZ  LIME  SHAPES 
MICHAEL  r.  PCTTER 

REAL  FUNCT  ION  TE  MPK  *8  ( W A VO  , S TRO  , ALF  AO  , ENGO  , OE  MOM,  TPWR , Ti  ,RTO, 

1 CTRFRO) 

IMPLICIT  REAL  * 81  A-H  , C'-Z) 

REAL  WAVO.STEO, AL^AO, ENGO, CTRFRO 

DATA  3LTZ , »LK,S PPL/ l. 33062 200F-23.6 .626 1960C E-34, 2 .99792500E ♦ 10/ 
EHCKT=(  ENG.3*S°0L*DLK *QtNOM ) 

STNGTH=STPT*TPWR*3EXP ( EHCKTM 1 .0-T l ) ) 

ALFA  = ALFA  0*RT0 
0 I F=CTR  FRQ—WA70 

SH AP c - ALFA/I (ALF A* ALFA)  + (CIF*CIF)> 

TE  MPK=SHA  PE*  STNG  TM 

RETURN 

ENO 


C 

C 

C 

C 

C 


130 


r 


00  01 
0002 
0003 
00  04 

00  )5 

00  06 

0007 

0008 

0009 

0010 
3011 


E OP  TRAN  IV  S LEVEL  21 


C 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE 
MICHAEL  T. 


MAIN 


PART 

POTTER 


DATE  = 77108 


11/01/56 


THIS  SUBROUTINE  COMPUTES 
VALUES  FOR  EACH  *'0L  EC  ULF 
LEAST  SQUARES  ?NC  DECREE 
TEMPERA  TURE . 


VIBPAT  IC.MAL  PARTITION  FUNCTION 
TYPE.  ARRAYS  4,3,C  ARE  DATA  FROM  A 
POLYNOMIAL.  DEPENDENT  VARIABLE  IS 


10 


SUBROUTINE  PART  (T.QRTO) 

IMPLICIT  3E4L*8( A-H.C-ZJ 
DIMENSION  A I 7) , P( 7) t C(  7) , QRTOI  71 

DATA  A/ 93 3 96 03? E-0 7, -.25483461 R-0 5 , - .2 167 3 75  IE-05 , 29 10 13726-05 

1..  00000 00  0,-. 60C6S938E-06,-. 98896 C39E-07 / 

DAT  A 3/ .445331  93E-94,  .516514206-03,  . 67 73 7 503E-03 , . 42 82 3000E-03 , 

1. 00000000.  .22766513E-03 ,.449331036-04/ 

DATA  C/  .99  504  4 24,  . 1C7C76  13E  01,.  5 65  6944  8,.  U2  01483E01  ,.  I 000000  0 E01  , 

1 .93548537 . . 99504424/ 

00  10  11=1,7 

QRTOI  T1  ) = A(  I1)«T*T  ♦ Bdlin  + CIIl) 

CONTINUE 

RETURN 

END 


SPACE=(80,  I 200,  ICC)  .RISE)  ,DCB=BLKSI  ZE  =80 


00100014 


XX 

//PQPT.SYSIN  OD  * 


IFF  2361 

ALLOC 

IEF237I 

449 

I EF  2 37 1 

432 

!EF  23  71 

236 

IEF237  I 

411 

I E F 1 42  I 

- STE 

I F p 28  5 I 

SYS 

I E F 2 35  I 

VOL 

UT'PllI 

STF  P 

JT :°l 3 I 

STFP 

STEPNAME 


PASSEO 


0094K 


> 


, FOR  J0.3N7838  FORT 
ALLOCATED  TO  SYSPRINT 
ALLOCATED  TO  SYSPIJNCH 
ALLOCATED  TO  SYSLIN 
ALLOCATED  TO  SYS  IN 
» WAS  EXECUTED  - CCND  COOP  0000 
77 ID8. T0Z5922.RVDC0. JOB N78 38 . LOAD SE T 
SER  NOS  = UT  F°  33 . 

FORT  EXECUTION  TIME  = 012.92  SEC. 

PORT  RFG1CN  = 0192K  MAX  CC°  E USED  « 

XXI.KF0  = XCC  PG  M*  | F W L » R E G i CN*96  K , P ARM*  (X  ,L  ET  . L 1ST ), COND=  1 V.LT.FOR’ 
DSNAME*SYS  l.FD? TL l 8 ,01  SP'SHR 

0SN=ijTEP.F3RTL  le,  DISP'SHR  FORTGCLG 

r>SN*,JTcP.CALCOMP.LGAOfCISP  = SHR  FORTGCLG 

SPACE* ( CYL , I 5, 5, 1 ) ) 

DSNAME'CGOSET I MA  I N ) , 0 IS  P=  I NEW  , P ASS  ) , UN  IT=SYSDA, 

SPACE' ( 1024,  < 20,10,1)  .RISE)  ,DC3  = «>LKS  I ZE  = 1024 
SYSOUT=A 

UNI T=SYSOA,SPAC t= (1D24, ( 1 00 , 1 0 ) , RL S E 1 , DC3= BL K S IZ E* 1 024, 
DSNAMt  = CSYSUT  1 

DSN  AME'CLOAOSET , 0 I SP= ( OL D, DEL ET E ) 

DONA  ME  = SY  SI N 

ALLOC.  FOR  J 0 BN  7 8 3 8 LKFD  STEPNAME 


XXSYSL IB 

DD 

XX 

OD 

XX 

DD 

D.SYSL  3 JO 

OD  1 

X/SYSLMOD 

DD 

XX 

XXSYSpp INT 

DO 

XXSYSUT1 

DD 

XX 

XX  SY  SL  I N 

DD 

XX 

DD 

00120000 

00140000 

00140020 

00140040 


XOO 160000 
00180014 
00200000 
X00210018 
00220018 
00240000 
00260000 


IEF236I 
ITF  ?37 I 
IE“  2371 
I E F 2 3 7 I 
IFCT37I 
IFF237I 
I E c 23  7 I 
IEP  2371 


131 

136 

136 

237 

44D 

232 

236 


ALLOCATED 
ALLOCA  TED 
ALLOCATED 
ALLOCATED 
ALLOCA  T=D 
ALL  C C AT  E D 
ALLOCA  TFO 


TO 

TO 

TO 

TO 

TO 

TO 

TO 


SYSL  13 


SYSLMOD 
SYSPR  INT 
SYSUT1 
SYSLIN 


i R = >42  I 

I E P ? 3 5 I 

- ct;i  WAS  EXECUTED  - COND  CODE  0300 

SYS1 . C0R  T L I 8 

KEPT 

I E p 285 I 

I EF??5  I 

VDL  SER  NO  S=  LNK213. 

UTEP  .FDRTL I 3 

KEPT 

I E F 7 35 1 
I-F235I 

VOL  SER  NOS  - UTEPll. 

UTEP.CALCOMP.LOAD 

KEPT 

IEF285I 

I Ec  23  51 

VOL  SER  NOS=  UTEPll. 

SYS  771 08. TD65922.R VOOO. JGBN7838  .GOSET 

PASSED 

I FF235  l 

I E F 28  5 I 

VOL  SER  NO  S=  UTEP  3 2. 

STFS77108  .T0b592  2 .RVDOO  . JOBN7838  . SYSUT  l 

DELETED 

IFP2  351 

I Fp285  T 

VOL  SFR  NO  S = OTESTl. 

SYS  77  108  .TD65922.RV330. JOBN 76 38. LOAD SE T 

DELETED 

I Fc  2351 

VOL  SER  NOS  = UTEP33. 

UTEP11I 
UT  = "13  I 


STEP  LKED  EXECUTION  TIME  = 002.04  SEC.  . 

STEP  LKEO  RFGICN  ' 3102K  MAX  CORE  USED  * 0132K 

XXGO  EXEC  PGM=*.  LKEO. SYSLMOD, COND'I 14 ,LT .FORT) , 14, LT.LKED) J 
XXFT  D5FD0 l DO  DONAMF=SYS IN 
//G0.‘ T06CC01  DD  SYS  OUT  = I A, ,6410) 

X/PT06F001  DD  SYSOUT=A 
XXFTJ7F0D1  DD  SYS(3UT  = B 
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